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ABSTRACT 


A computer program, surf3d, that uses the 3D finite element method to calculate the 
stress-intensity factors for surface, corner and embedded cracks in finite-thickness plates 
with and without circular holes, was developed. The cracks are assumed to be either 
elliptic or part-elliptic in shape. The computer program uses eight-noded hexahedral 
elements to model the solid. The program uses a skyline storage and solver. The stress- 
intensity factors are evaluated using the force method, the crack-opening displacement 
method and the 3D virtual crack closure methods. 

In the manual the input to and the output of the surf3d program are described. This 
manual also demonstrates the use of the program and describes the calculation of the 
stress-intensity factors. Several examples with sample data files are included with the 
manual. To facilitate modeling of the user’s crack configuration and loading, a companion 
program (a preprocesssor program) that generates the data for the surfSd called gensurf 
was also developed. The gensurf program is a three dimensional mesh generator program 
that requires minimal input and that builds a complete data file for surfSd. The program 
surfSd is operational on Unix machines such as CRAY Y-MP, CRAY-2, and Convex C-220. 
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INTRODUCTION 


Stress-intensity factors are fundamental quantities .used to predict fatigue crack prop- 
agation rates and crack growth profiles. Surface and corner cracks usually initiate at 
imperfections and voids in metallic structures. These cracks usually grow into near or 
part-elliptical shapes [ij. Therefore stress-intensity factors for elliptical cracks are needed. 
A computer program, surfSd, that uses 3D finite elements was developed to calculate the 
stress-intensity factors for surface and corner cracks in finite plates and in plates with 
circular holes. This program was used extensively by the authors [2-15] and the results 
were compared to those obtained by others and other methods. 

The purpose of this manual is to document this program, to describe the input to 
and the output of the program, and to demonstrate the program. Several examples 
are presented and several sample data files are included with this manual. To model 
crack configurations and loadings, a companion program (a preprocesssor program) that 
generates the data for the surfSd called gensurf [16] was also developed. The gensurf is 
a three dimensional mesh generator program that requires minimal input and that builds 
the complete data file for surfSd. 

First, the crack configurations and loading that can be analyzed with surf3d are dis- 
cussed. Next the program specifications and organization is presented. The procedure for 
the development of the models is explained. Then models for elliptic cracks are presented. 
Finally, several example problems and their output are presented. Appendix A lists names 
and functions of subroutines and major program variables. Appendix B describes the pro- 
cedures to compile and execute surfSd both interactively and in the batch mode on UNIX 
supercomputers such as CRAY Y-MP. 

CRACK CONFIGURATIONS AND LOADING 

Several crack configurations can be analyzed with surf3d. A cross-section through 
the crack plane of each configuration is shown in Figure 1. The configurations are (see 
Figure 1) 


(a) Surface crack 

(b) Embedded crack 

(c) Corner crack 

(d) Corner crack from a circular hole 

(e) Surface crack at a semicircular notch 

(f) Surface crack from a circular hole 

The cracks are assumed to be elliptic, semi-elliptic or quarter elliptic in shape and 
are defined by the semi-major axis, c, and the semi-minor axis, a. Any point on the crack 
front is defined by c,a, and </>, the parametric angle of the ellipse (see Figure 2). 

The first three cases (a), (b), and (c), can be analyzed by imposing appropriate 
boundary conditions on the model shown in Fig. 2(b). Similarly, the next three cases, 
(d), (e), and (f) can be analyzed by imposing appropriate boundary conditions on the 
model shown in Figure 2(c). The boundary conditions for all the six cases are described 
below. 
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For all the six cases, v = 0 is prescribed for all nodes on the uncracked portion (shaded 
portion in Figure 2), including the nodes on the crack front, of the y = 0 plane. 

(a) Surface crack in a plate: 

u = 0 for all nodes on the x = 0 plane 
w = 0 for the node at ( VF,h,0) 

(b) Embedded crack in a plate : 

u = 0 for all nodes on the x = 0 plane 
w = 0 for all nodes on the z — 0 plane 

(c) Corner crack in a plate : 

u = 0 for nodes at (0,h,0) and (0,/i, — t) 
w = 0 for the node at (W, h,0) 

(d) Corner crack from a circular hole : 

u = 0 for all nodes on the x — -R plane 
w = 0 for the node at (VF,h,0) 

(e) Surface crack from a semicircular hole : 

u = 0 for nodes at ( — R,h,0) and ( — R,h,—t) 
w = 0 for all nodes on the z = 0 plane 

(f) Surface crack from a circular hole : 

u = 0 for all nodes on the x = -R plane 
w = 0 for all nodes on the z — 0 plane 

Loading 

Several types of loading conditions can also be imposed on the cracked configurations. 
Four types of loading are commonly used: 

(a) Remote tensile loading 

(b) Remote bending loading about the x-axis 

(c) Remote bending loading about the z-axis 

(d) Uniform crack face pressure loading 

The three loading conditions (a), (b), and (c) are illustrated in Figure 3. These 
loadings are usually applied on the y = h plane as shown in Figure 3. The uniform 
pressure loading condition on the crack face is usually described by cr y = -1 at all the 
nodes on the crack face. 

The magnitudes of the tractions at the four nodes of a loaded face of an element are 
taken from the global vector of the nodal tractions. surf3d assumes that the tractions 
at the nodes are input in the global Cartesian coordinate directions. Using these nodal 
tractions, the tractions at any point on the loaded face of an element are obtained by 
interpolating with the element shape functions. From these tractions the consistent loads 
at the nodes are calculated following the standard FE principles [17]. 

Sometimes the user may require stress-intensity factors for loadings applied simulta- 
neously on several faces of the model. In these situations surf3d may be executed with 
(a) several loadings in a single run or (b) once for each loading condition. In the first 
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case, case (a), however, the boundary conditions should be the same for all the loading 
conditions and such that the loading produces only mode-I type deformations. 

In addition to traction type loadings, surfSd accepts displacement type loadings. The 
same restrictions mentioned in the previous paragraph for combined loading conditions 
apply. 


FINITE ELEMENT ANALYSIS 

The cracked solid is modeled with 8-noded hexahedra (Hex-8), isoparametric ele- 
ments. These isoparametric elements are serendipity elements based on a linear dis- 
placement field [17]. Near the crack front, singularity elements that are in the shape of 
pentahedrons are used. The formulation and the details of these elements can be found 
in references 2 through 4. 

A typical model of a quarter of a plate with a surface crack (with a/c = 1 and 
a/t = 0.5) is shown in Figure 4. At any station on the crack front, eight singularity 
elements are used as shown in Figure 4. Around the crack front these singularity elements 
form a torus. The rest of the solid is modeled with the Hex-8 elements. Models such as 
this were developed using the mesh generator program gensurf ; the details of the program 
are explained in the gensurf user’s manual [16]. 

Figure 5 shows the Hex-8 element and the pentahedral singularity element and how 
the elements are defined by their nodal connectivity. This figure also shows several ways 
to define the Hex-8 and singularity elements. The \/T terms in the shape functions require 
that T} is always along the crack front in the singularity element. 

Improper definition of the element connectivity leads to the calculation of zero or 
negative volumes. For example, two inconsistently defined elements are shown in Figure 
6. When an element’s volume is computed as zero or negative, surfSd aborts execution 
and prints the element number, nodal connectivity information, coordinates of each of the 
nodes of the element, and the value of the diagonal terms of the element stiffness matrix. 
The user can then use this information to check his input for improper definition of the 
elements. 


STRESS-INTENSITY FACTOR CALCULATIONS 

The stress-intensity factors are calculated at each station on the crack front using 
three methods. 

1. Force Method 

2. Crack Opening Displacement (COD) Method 

3. Three-Dimensional Virtual Crack Closure Technique (3D-VCCT) 

The first two methods are used when singularity elements are used at the crack front, 
while the 3D-VCCT method is used when the solid is modeled completely with nonsingular 
elements. 
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Force Method 


The force method was developed in references 2 and 3 and is briefly explained below. 
The method assumes that the 2 D state stress is valid within every infinitesimal portion 
of the crack front, so that the stress normal to the crack plane, <r y , can be written as 


a v - -7=== + A\ +0(r 1 / 2 ) 
V 2xr 


( 1 ) 


where Ki is the stress-intensity factor and Ai is a constant. Note that the distance r in 
Eq.(l) is measured normal to the crack front. The total force normal to crack plane and 
in a region bounded by z\ < z < z 2 and 0 < r < f£> is 


K 



o y dr] dz 


F y = -75 L2y/ro{z2 - zi) + AiTd(z2 ~ *i) + ■ • • 

V^Jr 

li 

— — j=2y/ro{zi - z{) + Air D (z 2 - z\) 
v2jt 


( 2 ) 


The forces F y in the region zi < z < z 2 and 0 < r < r D are known from the 
finite element analysis and can be used to evaluate the unknowns K\ and A\ in Eq. ( 2 ). 
However, Eq. (2) is valid only for small distances r from the crack front. A procedure for 
evaluating the if- values using the forces F y at various distances r is explained below. 

Figure 7 shows a portion of the crack plane with two consecutive layers i and (i + 1 ) 
in the model. The nodes j,k,...p and the nodes b,d,e,...h define the model ahead of the 
crack front in the i th layer. The FE solution calculates the forces F y at all these nodes as 
shown in Figure 8. In this figure, Fk t is the force in the y-direction computed at node 
k due to element I and F k y is the force in the y-direction computed at node k due to 
element J . The total force in the y-direction at node k in the i'Mayer is 


F k = F ki + F k2 ( 3 ) 

Thus the total force F y for r = r 2 in the i th layer (see Figure 8) is 

(F s ) r=ra = Fj + F k + F h +F b + F d + F ei ( 4 ) 

Similarly, the total force F y for r = r3 in the i th layer is 
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( 5 ) 


(•F v ) r _ rj - Fj + F k + Fi + F mi + Ft -f F 4 + F e + F ft 


Thus total F y forces are computed for five values of r and are used in Eq. (2) to obtain 
the following set of equations. 


IS 

(Fy)r=n = Fy , = —^=stilyfr\ + Ajritj 
{Fy)r=r 2 = F v , = —f=:ti2y/r2 + j 4 jrjtj 

{Fy)r=rs ~ F Vs — -^==ti2y/rs + i4ir s <i 


( 6 ) 


where i* is the thickness of the t th layer. The unknowns in Eq. (6) are Kj and A\. The five 
equations in Eq. (6) are used in a least square procedure to evaluate Ki and A\. Through 
numerical experimentation it was found that consistent results are obtained when 

(a) Five force equations are used (as in Eq. (6)), and 

(b) the maximum value of r(r 5 ) is less than or equal to (a/10), where a is the depth 
of elliptical crack (semi-minor axis). 

Equation (5) was modified slightly to calculate Ki at the two nodes (nodes j and b 
in Figures 7 and 8.) that define the ends of the crack in the i th layer. The total force on 
the ‘bottom’ * half of the * t,l layer, for example, for r = is 


(^»)r=r 2 - Fj + Fk + F t , = (F V2 ) ko< (7) 

Substituting Eq. (7) in Eq. (2) one obtains 

(K,L, = + M, )*,.■,(!) (8) 

* Note that ‘bottom* and ‘top* are relative and are used here for convenience in presentation. Each 
layer has a top and bottom. For example, node j in figure 7 is at the bottom of layer % and is also at the 
top of layer t - 1. Similarly, node b is at the top of layer land at the bottom of layer i + 1. Also note that 
the nodes j and b in figure 7 are at stations i and t + 1, respectively. 
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Note that this force is the force calculated using <r„ only in the bottom half of the i th 
layer. Similarly one can write for the top part of the * ,A layer as 

(F ” ) “' = ^ £ 2 2 ^ +(/l, W r2( l ) <»> 

where 


(■^y^top = Fb + Fd + F ei (10) 

Again the least square procedure is used independently for both top and bottom of the 
i th layer individually to calculate the stress-intensity factors at both sides of this layer. 

This procedure is repeated for all layers. The stress-intensity factors are then calcu- 
lated as 

( A/ ) statio, ‘ (') = 2 ^^ to P^layer ( i-1 ) bot ^layer (i) ^ 

for all stations except the first and the last. (Note that node j in Figure 7 defines station 
», node 6 defines station (i -f 1) , and so on.) For the first station, K[ is calculated from 

( A ') station (1) = i( K *) bot^ ayer ^ 

and for the last station, 

( A '') station (Nlaycr+1) = (Mayer; 

COD (Crack-Opening Displacement) Method 

In the COD method, the crack opening displacement at the nodes just behind the 
crack front are used to calculate the stress- intensity factor from 

COD = 2v = 2^^4(l-v 2 ) + 0{ r 3 ' 2 ) + ... (11) 

where G is the shear modulus and v is the Poisson’s ratio of the material. 
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The crack opening displacements are known at all nodes behind the cr ack front. 
Therefore, as in the force method, equation (11) is evaluated at several values of r as 
follows 


2Gv 


K J + A T + • • • 
Ki + A'r 


( 12 ) 


where Kj and A are unknown constants. For various values of r the nodal opening 
displacement v is known; hence, Eq. (12) can be written as 


£ ' fo)r=n 

VnT 

c • (i>) r=ra 

y/r* 


= K t + A'ti 


= Ki + A'r 2 


(13) 


£ • (tt)r=r» 

y/^S 


Kj + A r$ 


where c is a constant equal to 2 G y/2x / [4(1 — i/ 2 )]. 

Again the least square procedure is used to evaluate the constants Kj and A' in 
Eq. (13). As in the force method, five values of r are used in Eq. 13 and the nodes in the 
region 0 < r < (a/10) are used. The nodes that are used in the COD method k’ ,/',•••, p 
for the i th layer are shown in Figure 8. 

3D VCCT 

When the model does not contain singularity elements, the stress- intensity factors can 
be calculated from the strain energy release rates, Gj, assuming a state of plane strain 
[18,19]. 

For the i th layer, reference 19 proposed that 

(Gl)i = ~2Ka^ FiVk> + FhV ^ ( 14 ) 


where A,- is the average of the radial distance between the set of nodes j and k and nodes 
b and d, and t t is the thickness of the i th layer . The nodes k' and d' are the COD-nodes 
on the crack front (see Figure 8). Note that the forces Fj and Ft, are the force calculated 
using the elements in the i th layer alone. The stress-intensity factor is then calculated 
assuming plane strain conditions as 
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(15) 


Ki, = 


I EGj , 

(1 - j/ 2 ) 


The stress-intensity factor calculated from Eq. (15) is assumed to be the value at the 
center of the i th layer. 

A slightly different approach, suggested in reference 18, is used in surfSd and is 
described below. The node j on the crack front (see Figure 8) belongs to both layers 
(t- 1) and (»'). The strain energy release rate at node j (i.e., station i) in Figure 8 is given 
by 




(16) 


where 


= ) layer (^i)layer^ 

A = radial distance between nodes j and k 

= (also radial distance between nodes j and k ) 


(17) 


In Eq. (16) t Fj is the total force in the y-direction at node j , f a v is the average 
thickness of layers (i - 1) and (i) computed as in the force method. 

Eq. (16) is used for every station along the crack front. The stress-intensity factors 
are evaluated as before assuming plane strain as 



(1 - u2 ) 


or assuming plane stress as 


(A'/),- = ft • (<?,),- 


(18) 


( 19 ) 


Input Data for Stress-Intensity Factors 
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Figures 7 and 8 show the i ih and (* + l) th layers and the variables that are used 
to calculate the stress- intensity factors. Figure 7 shows the modeling with 4 singularity 
elements and shows the nodes, MNODE, and the elements, KELEM, used by the three 
methods for the i th layer. Note that in this figure, for clarity, the modeling detail behind 
crack front and above the crack plane are not shown. Figure 8 shows the forces that are 
used in the force method and the COD nodes that are used in the COD and the VCCT 
methods. In this figure, the subscript IR in the arrays FCENT and FTIP denotes the 
IR th loading condition (right hand side). 

Figure 9 shows the details on the z = 0 plane for a typical surface crack model. Here 
eight singularity elements are used. As explained in the gensurf user’s manual, first the 
z - o plane with a width of t units and height of h units is modeled. This is termed as 
the base model. Figure 9 shows only a part of the base model for a/t = 0.8 and a/c = 1. 
In this base model, 151 nodes and 128 elements were used. This figure also defines the 
input data for variables that are needed in the stress-intensity factor computations. Figure 
10 shows details on the crack plane (y = 0 plane) and very close to the crack front for 
a 4-layer model. This figure also shows the complete data that is needed to define all 
the variables that are necessary in the stress-intensity factor computations. The tedious 
preparation of the input data can be avoided by using a mesh generator like gensurf. 

The required input for calculating the stress-intensity factors using any of the three 
methods is described in the following section. By referring back to Figs. 7-9, the user’s 
understanding of input parameters will be greatly enhanced. 
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INPUT TO surfSd 


The required input data is described in this section. The data is created on a file, 
DFN, and is given an input to surf 3d. A complete data file for surf 3d could be generated 
using a mesh generator program like gensurf. Several sample input data files are attached. 
This section and the examples section together with Figures 6-8 will completely explain 
the input. (In this section, for convenience in presentation, the phrases cards, card sets, 
data sets, and lines are used interchangeably.) 


Card 

set 

Format 

Variable 

Description 

1 

20A4 

TITLE 

Title of the job 

2 

A6 

POUT 

Output Option. Input SHORT for short 
output and XLONG for long output option. 

3 +t 

f denotes free format 

EMOD,NU 

Young’s modulus and Poisson’s ratio 

4 

♦ 

NPOIN, NELEM 

Number of nodes and elements in the model 

5 

♦ 

I, X(I,1), X(I,2), 

Node number, z— , y— and z— coordinates of 


the model 


* X(I,3) 

Warning: surfSd assumes that Node 1 is at (c,0,0), i.e., at the end of the major 
axis of the elliptic or part-elliptic crack. 


i, 

Element number, nodal connectivity and the 

N0D(I,1), 

index of each element (=1 for singularity 

N0D(I,2), 

elements and =0 for Hex-8 elements). 

N0D(I,8), 


NINDX(I) 



Repeat this line NELEM times until all the elements are specified. 


11 







7 


NP,IU,IV,IW 


Boundary Conditions - 
Node number, 

Restraint code for the u— displacement 
Restraint code for the v- displacement 
Restraint code for the tv — displacement 
Restraint code: 

= 0 is free 
= 1 is fixed. 


Repeat this line until all boundary conditions are prescribed. 
Terminate this set of lines with four zeros. 


8 * 

NCASE 

Number of loading conditions (number of 
right hand sides) 

9 A6 

CTYPE 

Loading type. 

Input REMOTE for remote loading. Input 
CFACE for crack-face pressure loading. 

10 * 

NF, NL, NI, 
LIND, I FACE 

Pressure loading definition - 
Number of first element, 

Number of last element, 

Increment, 

Load index-input Unity, i.e LIND=1 
Number of 

face where pressure loading is prescribed. 

Repeat until all loaded elements are defined. 
Terminate this set of cards with five zeros. 

11 * 

IP, PX, PY, PZ 



Traction Definitions - 
Node number 

Magnitude of traction in the global x — direction at node IP 
Magnitude of traction in the global y - direction at node IP 
Magnitude of traction in the global 2 -direction at node IP 


Repeat this line until all the tractions are specified for current loading condition. 
Terminate each loading condition with zero for the integer and zeros for all three tractions. 
Repeat Item 11 until all traction loading conditions (NCASE) have been defined. 


12 



12 * NP,IU,IV,IW, 

U t v, W Nodal Displacements - 

Node Number, 

Restraint code for the u-displacement, 
Restraint code for the v-displacement, 
Restraint code for the w-displacement, 
Prescribed U-displacement, 

Prescribed V-displacement, 

Prescribed W-displacement 
Restraint code: 

= 0 is free 
=• 1 is fixed. 


Repeat this line until all displacements are prescribed. 
Terminate this set of lines with zeros for all four integers. 


13 * IRENUM 

Renumbering option: 

= 0 if no renumbering given; 

= 1 if renumbering scheme is provided. 


* JNEW(NPOIN) Renumbering scheme: Needed only if 

IRENUM is not equal to zero. 

Use as many lines as needed to read NPOIN integers. 


15 * NLOAD 


Number of concentrated loads 


lfi * NA(i) 
NA(2) 


Degree of freedom at which the concentrated 
loads are applied. 


NA(NLOAD) 


Use as many lines as needed to read NLOAD integers. 
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XY j 1 ) Magnitude of concentrated loads correspond- 

XY ( 2 ) * n 8 to the degree of freedom directions defined 

in NA(I). 

XY(NLOAD) 


Use as many lines as needed to read NLOAD values. 


18 * NSINGU, NLAYER 

Number of singularity elements in 
each layer and number of layers 
in the model. 


Note that the number of singularity elements are assumed to be the same in each layer of the model. 


19 * IC0D(I,1), 

IC0D(I,2), COD Node Definitions - 

IC0D(I,3), Node numbers of the nodes that are 

1C0D(I,4), used in the COD method (see Fie 81 

IC0D(I,5) ' 


Read (NLAYER+1) lines to define all the COD nodes along the complete crack front. 


MTIP(I,1), Node numbers defining the crack front for 

MTIP(I,2) each layer of the model with singularity 

elements (see Fig. 9). 


NTIP(1>1), Element numbers for pentahedral (singular- 

NTIP(I,2), ity) elements around the crack front for each 

layer of the model (see Fig. 9). 

NTIP(I, NSINGU) 


22 * KELEM(I,1), 

KELEM(I,2), 

KELEM(I,5) 


Element numbers of elements on the crack 
plane and ahead of the crack front for each 
layer of the model (see Figs. 7 and 9). 
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23 * MN0DE(I,1)> 

MN0DE(I,2), 

MNODE(I,10) 


Node numbers of nodes on the crack plane 
and ahead of the crack front for each layer of 
the model (see Figs. 7 and 9). 


HT, 

WIDTH, 

AOT, 

RPT, 

AOC 


Height of the model, h 
Width of the model, W 
( a/t ) ratio 

(i2 + t) value for a plate with a hole 
(a/c) ratio. 


As previously mentioned, Node 1 is assumed to be at (c,0,0). Thus the value of c is automatically 
defined. Using this value of c and the ratios AOC, (a/c), and AOT, {a/t), the values of a and t 
are computed by the program. 
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OUTPUT OF surf3d 


The output of sur/Sd is described in this section. Two output options lone and 
short, are available. If the user exercises the long output option (XLONG on the second 
card in columns 1-5), the following items will be printed. If the user specifies the short 
output option (SHORT on the second card), the items listed below with an asterisk will 
be omitted. Only a general description of major output sections is presented. Example 
output files are discussed and presented in the next section. 


Title 

Output Specification 
Description of the model 
Nodes and coordinates 

Element connectivity and index of the element. 

Boundary conditions 

Pressure or traction loading - elements, faces, nodes & traction magnitudes. 

Prescribed displacements 

Renumbering option and renumbering scheme 

Data for calculating stress-intensity factors 

Volume of the solid modeled 

For each loading condition: 

(a) Sum of x-forces before boundary conditions 

(b) Sum of y — forces before boundary conditions 

(c) Sum of z — forces before boundary conditions 


Projected surface areas: x-, y- t and z- components 

Nodal displacements at each node in the model for the first loading condition. 

•Nodal displacements at each node in the model for each loading condition. 
Nodal forces at each node in the model for each of the loading condition (If 
the absolute value of the forces F x> , F y> and F Zi at node t are less than 1CT 6 
then the forces are not listed.) 

•Average nodal stresses at each node for each loading condition. 

Equilibrium checks for each of the loading conditions 
(Ei F x , , Zi Fy, , E, F*. , * = 1, NPOIN in the model ) 

Sum of applied loads and surface area components for the current loading 
condition. ° 

Nominal stresses for the current loading condition. 

'Stress-intensity factor calculations:- 


COD method:- A'- value for each station along the crack front 
Force-method:- Apparent K - values at various distances from the crack 
front. 

These values are computed on ‘top’ and ‘bottom’ 
for each layer. 
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Summary of stress-intensity factors calculated using 

(a) the Force method 

(b) the COD method when singularity elements are present 

(c) the 3D-VCCT method, when singularity elements are not present 
in the model. K-values using both plane stress and plane strain 
assumptions are output with this method. 


Element Equilibrium: If each element satisfies equilibrium, £ F* - £ F y - 
y' p < j oJ5 — 6, then the program prints “All elements satisfy equilibrium . 
If some elements do not satisfy equilibrium, the program prints a warning 
message that N number of elements do not satisfy equilibrium and lists the 

element numbers. 


EXAMPLES 


In his section several examples illustrating the use of surfSd are presented Only 
part of the input data files are shown because the data files are very long. The cutout 
files are also very long even when the short output option is exercised, and, therefore only 
the pertinent parts of the output file are presented here. However, the complete data and 
output files are available on the disk accompanying this manual. 

In all examples the following assumptions are made: For remote loading the applied 
stress S is assumed to be unity. The maximum bending stresses (the outer surface fiber 
bending stresses) S bx and S bz are also assumed to be unity (see Figure 2). The crack 

depth, a, is assumed to be unity. The square of the elliptic integral of the second kind 
Q, is approximated in the program as ’ 


Q = 1.0 + 1.464(a/c) 165 , 
= 1.0 + 1.464(o/c) 1 ' 85 , 


for a/c < 1.0 
for a/c > 1.0 


All the stress intensity factors are normalized by S n ^5jjT where 


Sn = S, 

= s b „ 

= Sb, , 


for remote tension 

for bending about x — axis 

for bending about z-axis 


The loadings in these examples in this section are selected so that the nominal stress 
6 is unity. In general this will not be the case. Therefore, sur/Sd calculates the nominal 
stresses in the x- t y~, and z-directions for each of the loading conditions. The nominal 
stresses are calculated by evaluating the sum of the x-,y~, and z-components of the 
forces on the loaded faces of the model and dividing these components by the corresponding 
nonzero x-,y~, and z-components of the area of that face. The normalized stress 
intensity factors computed by the program can be then be divided by S where S is 
the correct value of the nominal stress. For bending loadings the nominal stress is zero. 

Therefore, for bending cases the outer-fiber stress is used as the nominal stress. This is 
assumed to be unity in the program. 

For loading in the form of the prescribed displacements, the user needs to calculate 
the nominal stress. This is because the information on the loaded face areas are not 
available to surf3d. Fortunately, this calculation is very simple to perform. The sum of 
and ‘-components of the forces on this face is available from the output of 
surfSd. These force components are divided by the corresponding nonzero components of 
the area to obtain the value of S. The normalized stress-intensity factors given by the 
program are then divided by S to obtain the correct normalized value. Note that this 

division needs to be performed by the use externally, i.e., after obtaining the output from 
surf 3d. 

The following table defines the input and output files used in each of the examples 
presented in this section. v 
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Input and Output Files in the Examples 


Example 

Number 

Input 

File 

Output 

File 

1. 

dexl 

outl2 


(Table It) 

(Table 2 f ) 

2. 


outnl2 
(Table 3) 

3(a). 

dexSa 

outr22 


(Table 4) 

(Table 5) 

(b). 

dexS b 

outc22 
(Table 6) 

4(a). 

dexja 

outr28 

(Table 7) 

(b). 

dex4b 

outc28 

(Table 8) 

5. 

dex5 

outcor28 
(Table 9) 

6. 

dex6 

outem28 

(Table 10) 

7. 

dexl 

occorlS 
(Table 11) 

8. 

dex8 

oscorl5 
(Table 12) 

9. 

dex9 

osmcor!5 
(Table 13) 

10. 

dexlO 

outdl2 


(Table 14) 

(Table 15) 

11. 

dexll 

outdx!2 

(Table 16; 


t Partial listing is shown in these tables 
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Traction- Type Loading 

Surface, Corner and Embedded Cracks in Finite Plates 

Example 1: Surface crack with a/c = 1 and a/t = 0.2 in a plate subjected 

to remote tension and bending about x— and z-axes. 

(W = 25; /i = 125;/ = 5.0) 

One quarter of the solid is modeled with 2161 nodes and 1664 elements. Eight layers 
are used to model the solid. In each layer 8 singularity elements are used at the crack 
front. The input data file is partially presented in Table 1. Part of the output file is shown 
in Table 2. The complete output file outl2 is available on the disk. The first part of the 
output is for remote tensile loading (loading number 1). The stress-intensity factors at 
each station along the crack front from </> = 0 to r/2 calculated by the force and COD 
methods are given. In this part of the output, first the absolute value of K and then the 
normalized value of K from both methods are presented. Next the normalized values of 
K by the force and COD methods are tabulated for each station along the crack front, 
i.e., for each value of <f>. The results for loadings 2 and 3, bending about x- and z-axes, 
respectively, are also listed in Table 2. 

Example 2: Same configuration as in example 1 but without singularity elements. 


This example uses non-singular elements through out the model. The program scans 
the NINDX array to determine if singularity elements are present in the model. When 
singularity elements are not present, force and COD methods are not used and instead the 
3D-VCCT method is used to calculate the stress- intensity factors. This example presents 
the results obtained using the same model as in example 1, and setting all the indices 
INDX, in the NINDX array to zero. When 3D-VCCT method is used, one has to assume 
either plane stress or plane strain. The program calculates the normalized values of the 
stress-intensity factors using both assumptions. The results for example 2 are presented 
in Table 3. Comparison of these results with those in Table 2 shows that the 3D-VCCT 
assuming plane strain agrees well with the force method results for most of the crack front 
for all three loading conditions. 

Example 3: Surface crack with a/c = 0.2 and a/t = 0.2 in a plate subjected to 

(a) remote uniform tensile loading and 

(b) uniform crack face pressure loading. 

As in the previous examples, one quarter of the solid is modeled with 2441 nodes 
and 1872 elements. Again eight layers are used and in each of these layers 8 singularity 
elements are used at the crack front. Table 4 presents a partial listing of the data file. 

Prom superposition principles, it can be shown that the stress-intensity factors for 
the two loading conditions are identical. This example demonstrates that surfSd nearly 
reproduces this result. Tables 5 and 6 present, respectively, the partial output files for 
remote tensile and crack face pressure loadings. Comparison of the results shows that the 
stress- in tensity factors calculated by both methods are nearly but not exactly identical. 
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The slight differences between the two sets of stress-intensity factors are d,.P to tk„ 
exact nature of the loading on the crack face. Note that on the c,„„ tL ... n ° n ’ 

semi-elliptic slit with a semi-major axis of Jc 1 ~ (Q OOla^ anH B ^ €re ls a 

The loading for the two cases are not ^ ^ ° °' 0010 - 

differences in the stress-intensity factors. **** “ d * hence ’ the 

Example 4: Deep surface crack in a plate with a/c = 0.2 and a/t = 0 8- 

(W = 50;/i = 125; < = 1.25;o = 1 ; c = 5 ) 

One quarter of the solid is modeled with 2464 nodes anW i««; „i . . , , 

& «o e f ridt froni 8 s “ y “ ■ - *«s 

(a) remote uniform tensile loading, and 

(b) uniform crack face pressure loading of magnitude S. 

The input files for these two cases dex4(a) and (b) are on the disk T a ki D 7 a 
8 present the output files for loadings (a) and (b) respectively Cnmn • 7 ^ 

obtained for this deep elliptic crack show that both loadings give nearly identic^ 7^ 
intensity factors along most of the crack front g g y d nt,caJ stress ' 


Example 6: Corner crack with a/c = 0.2 and a/t = 0.8* 

(W = 50;h = 125;/ = 1.25; a = l;c = 5) 

Consider a quarter-elliptical corner crack in a plate subjected to remote tensile load 

identSaU*^ °! “ m ° dcled with 2464 »odes and 1856 elements. This model is 

identical to that use< I in example 4 but the boundary conditions are changed to the corner 
crack boundary conditions. The output is presented in Table 9. 

Example 6: Embedded crack with a/c = 0.2 and a/t = 0.8 

(ir = 50;/i= 125;/= 1.25;a= l;c = 5) 

is kuSf!" U,C S f d “ inodeled with 2464 and 1856 elements. This model 
identical to those used in examples 4 and 5 but the boundary conditions are changed 

to embedded crack-boundary conditions. The output is presented in Table 10. 6 

Surface and Corner Cracks in a Plate with Holes 

In this subsection, three examples of surface and corner cracks at a semi-circular or a 

circular hole arc presented. In all these examples, the crack front is modeled with 8 layers 

with 8 singularity elements at the crack front in each layer. The hole Ts Z^ed (in the 

f" °D T’ See ( tlgnr ° 2 \ w,th 5 unef l ual thickness layers. The thickness of the five layers 

hoi /ql P , P are 5,1 °’ 25 ’ 35 ’ and 40 P ercent the radius of the 

hole. (See gensurf users manual [16] for details on hole modeling.) 


Example 7: 


Corner era. :ks at a circular hole with a/c = l, a/t = 0.5, and R/t = V 
\W — 25; li = 125; a — 1; c = 1 ; / = 2; i? = 2) 
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Because of symmetries, one quarter of the solid is modeled with 2863 nodes and 2260 
elements. Loading conditions of remote uniform tension arid bending about x-axis are 
considered. Table 11 presents a partial listing of the output. 

Example 8: Surface crack at a circular hole with a/c = \,a/t = 0.5, and R/t = 1; 

(W = 25; h = 125; a = c = 1; R = t - 2) 

The same model as in example 7 is used for this example. The boundary conditions 
on r = 0 plane are prescribed to rellect the symmetric boundary conditions, i.e., all nodes 
on z = 0 Plane are prescribed to have zero ui-displacements. Remote uniform tensile 
loading was prescribed. Table 12 presents a partial listing of the output. 


Example 0: 


Surface crack at a semi-circular hole with a/c - l,a/t - 0.5, and 
R/t = 1- (W = 25; ft = 125; o = c = 1;H = t = 2) 

The same model aL in examples 7 and 8 is used in this example. All the boundary 
conditions, except those on the x = -R plane see Fig 2(c)), used in example 8 are used. 
Only one node on the x = -R plane is prescribed to have a zero u- displacement Remote 
uniform tensile loading was applied to this model. Table 13 presents a partial bsting of 

the output. 

Prescribed Displacement Loadings 


Example 10: 


Surface Crack in a Finite Plate 

Surface crack with a/c = 1 and a/t = 0.2 in a plate subjected 
to remote prescribed displacements, v = 10 in. on the y - h face. 
(W = 25; /» = 125; t = 5.0) 


The configuration in this example is identical to example 1 Part of the input data file 
is presented in Table 14. Part of the output file is shown in Table 15. The complete ou pu 
fill outdl2 is available on the disk. The calculation of the nominal stress is ^strated 
in this example. The sum of the forces on the loaded face y - h are compute 
F =00 F = 29.99872 lbs, and F z = 0.0. The components of the area are A* _ 
0 0 A v = 25 V 0x5.0= 125.0 sq. in., and A z = 0.0. Therefore, the nominal stress n .5 - 
29 99872/125.0= 0.24 psi. All the output stress-intensity factors are divided by 0.24, th 
v!iufof he nominal stress 5. For example, the K/S^jQ at * = */2 using the force 
method will now be equal to 0.24519/0.24 = 1.0216. This value agrees extremely well 
with that for the traction-type loading in example l for the same configuration (using 
force method this value is 1.0221). This is expected because of the large plate used in 

both examples. 


Example 11: 


Surface crack with a/c = 1 and a/t = 0.2 in a plate subjected 
to remote prescribed displacements, u - -0.3 • 10 in. on the x 
face. 

( W = 25; h = 125; t = 5.0) 


The configuration in this example is identical to that used in example 1. The loading 
used! thTexample gives zero values for the stress-intensity factors. This example shows 
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that the stress-intensity factors calculated by surf3d will be nearly zero but not identically 
zero. Note that the run was aborted because of the attempted square root of small negative 
stress- in tensity factors. Part of the output file is shown in Table 16. The complete output 
file outdxl2 is available on the disk. 
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APPENDIX A 


PARAMETERS, SUBROUTINES, MAJOR PROGRAM 
VARIABLES, AND COMMON BLOCKS 


This appendix presents the names and functions of the subroutines and major pro- 
gram variables and the common blocks with their elements and the subroutines which use 
the common blocks. A flow chart of surf 3d is presented in Figure 11. 

A-ls PARAMETER STATEMENT VARIABLES 


NAME 

DEFINITION 

MAXBK 

Maximum dimension of the assembled stiffness matrix - BIGK 

MAXNOD 

Maximum number of nodes 

MAXEL 

Maximum number of elements 

MAXRHS 

Maximum number number of right hand sides (loading condi- 
tions) 

MAXBC 

Maximum number of boundary conditions 

MAXB 

Maximum bandwidth of equations 

NNODE 

Number of nodes on the Hex-8 element = 8 

NDOF= NFR.EE 

Number of degrees of freedom per node = 3 

NSIF 

Number of stations along the crack front where the stress- 
intensity factors are evaluated. 

MAXRUND 

Maximum dimension of IRUND 

MAXDIS 

Maximum number of degrees of freedom 
= MAXNOD*NFREE 

NSMK 

Dimension of the stiffness matrix of the Hex-8 element. 
= NNODE * NFREE =24 


NA ME 

i. ADJUST 


2. ASEMB 

3. ASTAR 


A-2: SUBROUTINES 
FUNCTION 

Reorders the nodal coordinate array, the element connectivity 
array and the boundary condition array according to the 
renumbering scheme provided by the user. If no renumbering is 
given, no reordering of these arrays is performed. 

Processes all elements, obtains element stiffness matrices, load 
vectors and assembles the global stiffness matrix. 

Obtains the transformation relationship between the generalized 
coordinates and the nodal coordinates. This routine is called 
only once. 
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4. BLOCK DATA 

5. BOUND 

6. CCLOCK 

7. ODER 

8. CORDIN 

9. DERIVE 

10. FORCES 

11. CDERV 

12. LITTLE 

13. LOAD 

14. MATIN V 

15. MATMUL 

16. MODULUS 

17. PARENT 

18. REND 

19. SETPP 


Contains the Gaussian coordinates and weights up to an 8-point 
integration rule. Also contains the parent coordinates 
of the Hex-8 elements. 

Prescribes the boundary conditions. 

Calculates the CPU and accumulated CPU times between 
successive calls. 

Calculates the Cartesian derivatives at the Gaussian points and 
forms the BJ matrix. 

Reads the coordinates of all the nodes in the model. 

Obtains the derivatives of the shape functions at any point 
in the Hex-8 and pentahedron singularity elements. 

Calculates the element forces and checks the element equilib- 
rium. Calculates the nodal stresses and forces and checks global 
equilibrium. 

Calculates the parent derivatives, i = 1,8, at each 

of the ( NGAUSS) Gaussian points, both for singularity and 
Hex 8 elements (If reduced integration is used (IRED=1), the 
parent derivatives are calculated at the center of the element). 

Calculates the numerically smallest node number on each 
element 

Calculates the consistent loads on each of the loaded elements. 

Obtains the inverse of a square matrix. 

Obtains the product of two matrices. 

Computes the modulus matrix of an isotopic material. 

Calculates the parent derivatives, = 1,8, at each of 

the (NGAUSS) 2 Gaussian points. 

Calculates the half-band width for each degree of freedom, 
sets up the row pointer array and calculates the total memory 
required to store the assembled stiffness matrix in profile form. 

Calculates the half- band widths required for each of the degree 
of freedom in the model. 
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20. SHAPE 

21. SFACTOR 

22. SMALL 

23. SMALLK 


24. SOLV 

25. STRAN 

26. STRESS 

27. SYMBAN 

28. TRANF 

29. TRANS 

30. ZEROLN 

31. ZEROLV 


Calculates the shape functions of the Hex-8 element. 

Calculates the stress-intensity factors using the force and the 
crack-opening displacement methods. 

Controlling subprogram which calls the stiffness matrix and 
load routines for the Hex 8 and singularity elements. 

Calculates the element stiffness matrix of the Hex 8 and 
singularity elements. 

Note: The element stiffness matrices are calculated using a 
procedure similar to that presented by Noor and Hartley [20] 
for vector computers like the STAR-100, CYBER 203, CYBER 
205. The vector version was very efficient and was later 
devectorized for UNIX machines. Numerical experimentation 
showed that the current UNIX version is as efficient as the 
CYBER 205 version. The source code in this subprogram does 
not follow conventional procedures used for evaluating element 
stiffness matrices. 

Calls the solver, SYMBAN, and prints the displacements of the 
all the nodes in the model. 

Calculates the stress- transformation matrices for the singularity 
elements. 

Calculates the stresses in each element by calculating the 
stresses at the 2x2x2 Gaussian points. 

Solves the system of equations using the Cholesky decomposi- 
tion. The left hand side matrix is arranged in profile form. This 
is an in-core solver. 

Calculates the singularity element transformation matrices. 
Obtains the transpose of a matrix. 

Zeros out an integer array. 

Zeros out a real variable array. 
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A-3: MAJOR PROGRAM VARIABLES 


The variables in various commons arc listed alphabetically by the common block 
name and then the other major variables are listed in alphabetical order. 

COMMON VARIABLE DEFINITION 


ASTIF BIGK(MAXBK) 
AVERAGE FS UMS(3,MA XRIIS) 

ASUMS(3) 

BANDW P(MAXI)IS) 

T(MAXB) 

BNOD X(MAXNOD,NFUEE) 

CLOCKS TSTART 

COD ICOD(NSIF,MAXKE) 

CLIST LIST(MAXBC) 

NA(NLOAD) 

X Y(NLOAD) 

NLOAD 

COMB NCASE 


Assembled stiffness matrix stored in 
profile form. 

Sum of the forces in the x-,y- and 
^-directions on the loaded faces of the 
model. 

Projected areas in the x-,y- and z- 
coordinate directions of the loaded 
fares, 

Work array used to store the reciprocal 
of the diagonal coefficients in the solver 
Work array needed to store the half- 
band widths in the solver 

■r — ,y— and z— coordinates 
of all the nodes in the model. 

X(1, 1), X(I,2), X(I,3) represent the 
x-,y- and 2 — coordinates of node I. 

Temporary storage for CPU Time 

Nodes that are used in the COD 
method to evaluate the stress-intensity 
factors at various stations along the 
crack front. 

Iloundary condition array. LIST(I) de- 
fines degrees of freedom (dof) I to be 
zero. 

Defines the dof where external loads are 
prescribed. 

Defines the magnitude of the external 
load corresponding to the dof defined in 
the the NA array. 

Number of external loads prescribed. 
Number of loading conditions. 
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DEPTH 

NSINGU 


DISP 


GAUSS 


INI) 


Thickness of the solid (t, see Figure 1) 
Number of singularity elements around 
the crack front 

R HOLE Radius of the circular hole 

DIS(MAXDIS,MAXRHS) Array that contains the loads corre- 

sponding to each of the dof before so- 
lution and the displacements after the 
solution for each loading condition. 

CORD (8,8) Gaussian coordinates 

WEIGHT(8,8) Gaussian weights. 

An N- point Gaussian coordinate and 
weight can be found in the Nth column 
of the arrays CORD and WEIGHT, re- 
spectively, (1 < N < 8). 


NINDX(MAXEL) 


Array used to index the elements 
NINDX(I)=0 denotes that the 

element I is a Hex-8 
element. 

N IN DX( 1)= 1 denotes that the 
element I is a 
pentahedron singularity 
element. 


LINDX(MAXEL,2) 


Array used to indicate elements with 
pressure loading. 

LINDX(I,1)=0 denotes that the 

element I does not 
have any pressure loading. 

L1NDX(I,1)=1 denotes that the 

element I is subjected 
to pressure loading. 


LINDX(I,2) defines the face on which 
the pressure loading is prescribed. The 
6 faces are defined using the parent co- 
ordinates £,tj,C 18 follows: 


£=o 

face number = 

1 

{=1 

face number = 

2 

Tf—0 

face number = 

3 

r/-l 

face number = 

4 

<=o 

face number = 

5 

C=i 

face number = 

6 
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INTGR 


INTNST 


MASTER 


LIND 

INDX 

IFACE 

NPOIN 

NBOUN 

NELEM 

NBAND 

NDIS 


KELEM(NSIF,MAXKE) 


MN0UE(NS1F,2*MAXKE) 


FCENT(NSIF, 

2*MAXKE,2, 

MAXRHS) 


NKO UNT(NSIF,2*MAXKE ) 


MTIP(NSIF,2) 


NTIP(NSIF,MSINGU) 


FTIP(NSIF,2,MAXRIIS) 
N LAYER 


LIND = LINDX (1,1) 

INDX = NINDX(I) 

1FACE= LINDX(I,2) 

Number of nodes in the model 
Number of boundary conditions 
Number of elements in the model 
Maximum half-bandwidth of the 
equations 

Number of dof in the model 
= NPOIN*NFREE 

Elements on the crack plane and ahead 
of the crack plane used in stress- 
intensity factor calculation. (See Figure 
7 for definitions). 

Nodes on the crack plane and ahead 
of the crack plane used in the stress- 
intensity factor calculation. (See Figure 
7 for definitions). 

Forces in the y-direction at the 
nodes on the crack plane and ahead of 
the crack used in the stress-intensity 
factor calculation. (See Figure 8 for def- 
inition of this array.) 

An integer array to keep the count of 
the occurrence of nodes in the MNODE 
array. 

Nodes on the crack front of each layer of 
the model. 

(See Figures 8 through 10 for definitions 
and examples). 

The elements closest to and around the 
crack front. These elements are used to 
evaluate the forces at the crack front 
nodes. 

(See Figures 8 through 10 for definitions 
and examples). 

Forces at the crack tip nodes. 

(See Figures 8 through 10 for definitions 
and examples). 

Number of layers (wedges) in the model. 
Stress-intensity factors are evaluated 
at (NLAYER-f 1) stations on the crack 
front. 

(See subprograms GDERV and PAR- 
ENT for details) 
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DMASTX(216,3), 

DMASTE(216,3), 

DMASTZ(216,3) 

Derivatives for a 

Hex-8 element, at the quadrature 
points for 2 or 3-point Gaussian 
quadrature. 

For a 2-point Gaussian, 8 derivatives, 
For a 3-point Gaussian, 27 derivatives. 


DMASSX(216,3), 

DMASSE(216,3), 

DMASSZ(216,3) 

Derivatives for 

pentahedron singularity elements at the 
(8 or 27) quadrature points with 2- or 
3-point Gaussian quadrature. 


DS1IX(216,3), 

DSHE(216,3), 

DSHZ(216,3) 

Derivatives ^ for Hex-8 

elements at the center of the element 
for reduced integration 


WT(27,3) 

Gaussian weights at the quadrature 
points for 2- or 3-point Gaussian 
quadrature in each direction. 

PDER 

DNX(216), 

DNY(216), 

DNZ(216) 

Parent derivatives 
of all shape functions (*=1,8) at all 
the quadrature points (8 or 27) for 2- 
and 3- point Gaussian quadrature 
(see subprogram GDERV for details) 

POINTER 

NDW(MAXDIS) 

IRPNT(MAXDIS) 

NDSTK(MAXNOD) 

Bandwidth for each degree of freedom 
Row pointer array for each degree of 
freedom 

NDSTK(I) gives the lowest node num- 
ber that is connected to node I. This 
array is used to compute the IRPNT 
array in subprogram REND. 

RENUM 

JOLD(MAXNOI)) 

JNEW(MAXNOD) 

Array which relates the old nodal num- 
bers to the new nodal numbers in the 
renumbered scheme. JOLD(IN) gives 
the old number of new node IN. This 
array is complementary to the array 
JNEW. 

Array which relates the new node 
numbers to the old node numbers. 
JNEW(IO) gives the new nodal num- 
ber of the old node 10. This array is 
complementary to the array JOLD. 
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STIF NGAUSS Order of Gaussian integration used to 

integrate the elements 

I RED Flag used to turn on and off reduced 

integration 

IRED =1 Reduced 

integration is ON 
= 0 Reduced 

integration is OFF 


DELVOL 

SUMM STRESS-INTENSITY FACTORS 

SIFF(NSIF) 

SIFD(NSIF) 

TYRE CTYTE 

ULOAD UL(NSMK,MAXRHS) 

ALOAD( MAXNOD, 
NFREE,MAXRIIS) 

TLOAD(NNODE, 

NFREE,MAXRHS) 


Volume of the element being processed 

Stress- intensity factors calculated by the 
force method. 

Stress-intensity factors calculated by the 
COD method assuming plane strain. 

Loading Type: 

REMOTE-for remote loading 
CFACE - for crack face pressure loading 

Consistent loads calculated by subpro- 
gram LOAD for each of the loading 
conditions 

Magnitudes of tractions at each 
of the nodes in the model 
for all loading conditions 
Magnitudes of tractions at each 
of the nodes of the element being 
processed for all loading conditions. 


MAJOR PROGRAM VARIABLES NOT IN ANY COMMON BLOCK(S) 

VARIABLE DEFINITION 

E L DIS( NSM K, MA XRUS) Nodal displacements of the element be- 

ing processed. 

F0R(NSMK,MAXR1IS) Nodal forces of the element being 

processed. 

FORCE(MAXNOD,3,MAXRHS) Nodal forces at each node of the model. 

The first subscript indicates the node 
number, the second subscript indicates 
the direction, the third subscript indi- 
cates the loading condition. 
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ND(MAXNOD) 


Array containing the number of ele- 
ments connected to a specific node. 


SP(MAXN0D,6,MAXRHS) Average nodal stresses at each node of 

the model. The first subscript indicates 
the node number. The second sub- 
script indicates one of the six stresses, 

* °y * °* * °xy » Oyx > &xx- The third sub- 
script indicates the loading condition. 

ND(MAXNOD) Array containing the number of ele- 

ments connected to a specific node. 

XE(NNODE,3) Nodal coordinates of the element being 

processed. 


34 


APPENDIX B 
EXECUTION OF surfSd 


This appendix describes the procedures to compile and execute surfSd both inter- 
actively and in the batch mode on a CRAY Y-MP computer. To distinguish among 
commands and responses in this appendix the following code is used. The typewriter 
font is used to denote a system response. (For example, flyer 25X denotes the 25th 
command to the Cray Y-MP named flyer.) The italic font is used to denote the user 
commands and the file names. The Roman font is used in the explanation of various 
commands and responses. 

B-l: Interactive Compilation and Execution of surfSd 


flyer 25X 

cfm surf 3d. f 

Invoke the FORTRAN compiler 
and create the object file, 
surfSd. o 


FF0001 CFT77 VERSION 4.0.3.13 

(392409) 04/30/91 13:15:30 

FF0002 COMPILE TIME 14.168 SECONDS 

FF0006 MAXIMUM FIELD LENGTH 409116 

DECIMAL WORDS 

FF0003 4963 SOURCE LINES 

FF0004 0 ERRORS 

FF0005 CODE: 9683 WORDS, DATA: 2338 
WORDS 

( For Convex computers use 

Here - cfc flag emulates the 


fc -cfc -72 -o s.e surf 3d. f ) 

Cray FORTRAN compiler and 
the -o s.e names the 
executable s.e. If -o s.e is 
omitted the executable name 
defaults to a. out 

flyer 26X 

scgldr -f indef -o s.e surfSd. o 

link and create the 
executable file, s.e. 

The core is set to 
indefinite values with the 
-/ indef flag. 



flyer 27X 

s.e < datl2 >k outl2k 

Execution with datl2 as data 
file and outlS as output file. 

The first ampersand directs that 
the execution to be carried out 
in the background. The second 
ampersand directs the errors 
to the output file, outl2. 

flyer 28% 

ps 

Determine process status 
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PID TTY TIME COMMAND 
41299 pOOl 0:00 ps 
41253 pOOl 0:00 quotation 
41292 pOOl 0:14 s.e 
41247 pOOl 0:00 csh 

flyar 29X pa Process status 

PID TTY TIME COMMAND 
41253 pOOl 0:00 quotamon 
41315 pOOl 0:00 ps 
41292 pOOl 0:30 s.a 
41247 pOOl 0:00 csh 

flyar 30* 

[1] Done s.e < datl2 >t outl2t Job completed. 

B-2: Batch Compilation and Execution of surfSd 


The following commands describe how to compile and execute surfSd in the batch 
mode. A batch file, surfer , is used to perform the majority of this task. Because of 
variations in the batch processes in different systems, the users are encouraged to consult 
their system administrator to determine specifics of their system batch processes. 

Here the source and the data files are assumed to be in a temporary (space) directory 
/tmp/raju. 

flyer 25* cat surfer Lists the script file, surfer. 

#QSUB -me 
#QSUB -IT 100 
#QSUB -IM 16mw 

jo- 

set verbose 

date 

ed /tmp/raju 

pwd 

eft 77 surf 3d. f 

segldr -f indef -o s.e surfSd.o 
./s.e < datl2 > outl2 

date 

ja -c sft 

flyer 26X qsub -eo -o slog surfer Submits the script surfer 

to the batch processor with a 
request to write the log, 
the errors and all the 
information to the file 

slog. 
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flyer 27X 


qstnt -a 


Tells the user the status 
of batch job. 

Lists the log file slog. 


flyer 28X cat slog 

The table B-l lists only the pertinent portions of the log file. 



Warning: no access to tty; thus no job control in this shell... 


Beginning Batch Execution 

18:41:21 

18:41:21 ♦ exit 

18:41:21 + /usr/spool/nqs/scripts/++02f+++++0+*+ 

date 

Sun Sep 29 18:41:21 EDT 1991 

cd /tmp/raju 

/tmp/raju 

pwd 

/tmp/raju 
cft77 surf 3d. f 


segldr -f indef -o s.e surf 3d. o 
./s.e < datl2 > out 12 
STOP (called by $MAIN ) 

CP: 29.470s, Wallclock: 29.823s, 24.7. 

HWM mem: 7730286, HWM stack: 310502, Stack overflows: 0 

ja -csft 

Job Accounting - Command Report 
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Command Started Elapsed User CPU Sys CPU I/O Wait I/O Wait 
Name At Seconds Seconds Seconds Sec Lck Sec Unlck SBU’s 


ja 18:41:21 0.0042 0.0004 0.0037 0.0000 0.0000 0.00 

date 18:41:21 0.0026 0.0005 0.0020 0.0000 0.0000 0.00 
pud 18:41:21 0.0049 0.0005 0.0043 0.0000 0.0000 0.00 
pwd 18:41:21 0.0048 0.0005 0.0042 0.0000 0.0000 0.00 
cft77 18:41:21 14.3477 14.1998 0.1284 0.0030 0.0147 7.16 
segldr 18:41:35 0.6375 0.5865 0.0503 0.0000 0.0000 0.32 
s.e 18:41:36 29.8828 29.4708 0.3972 0.0007 0.0112 14.93 
date 18:42:06 0.0026 0.0005 0.0020 0.0000 0.0000 0.00 
Job Accounting • Command Flow Report 


parent -> child . . . 


ja 

date 

pwd 

pwd 

cft77 

segldr 

s.e 

date 

Job Accounting - Summary Report 


Job Accounting File Name : /tmp/nqs.+++*+002f/. jacct5018 

Operating System : snl015 flyer 6.0 woo. 19 CRAY Y-HP 

User Name (ID) : raju (14006) 
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Group Name (ID) : ncedu (14000) 

Account Name (ID) : raju (14006) 

Job Name (ID) : surfer (5018) 

Report Starts : 09/29/91 18:41:21 

Report Ends : 09/29/91 18:42:06 

Elapsed Time : 45 Seconds 

User CPU Time : 44.2593 Seconds 

System CPU Time : 0.5922 Seconds 

1/0 Wait Time (Locked) : 0.0037 Seconds 

1/0 Wait Time (Unlocked) : 0.0259 Seconds 

CPU Time Memory Integral : 231.3400 Mword-seconds 

SDS Time Memory Integral : 0.0000 Mword-seconds 

1/0 Wait Time Memory Integral : 0.0070 Mword-seconds 

Data Transferred : 3.1413 MWords 

Maximum memory used : 7918592 Words 

Logical 1/0 Requests : 1796 

Physical 1/0 Requests : 755 

Number of Commands 8 
Billing Units : 22.4258 

logout 

18:42:06 + clear 
18:42:06 + set lo=logout 
18:42:06 ♦ set bye*logout 
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B-3: Lower to Upper Case Conversions 


After the input to surfSd is generated, change all the lower case characters in the file to 
the upper case characters using the script file called trans. The conversion is recommended 
because of string inputs. When alphanumeric strings with mixed upper and lower case 
characters are used, comparisons cause problems. To avoid these difficulties all upper case 
characters are recommended in the input files. The conversions file trans is provided with 
this manual. For example, to change all the characters in an input file called tinp to upper 
case characters, type 

trans tinp 


The system response will be 

remove tinp.n? 

Type y to remove temporary files. All the lower case characters in the file tinp are now 
changed to upper case characters. If the original file does not contain any lower case 
characters then the above command has no effect on the file. Therefore, the use of trans 
is recommended on the input file before the file is submitted for execution with surfSd. 

The script file trans contains the following statements. 

• 

tr a-z A-Z <$1> $1 .n 
cp $l.n $1 
rm $1 ,n 
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(a) Surface crack In a plate 



(c) Corner crack In a plate 



(e) Surface crack at a 
semi-circular notch In a plate 


T 

t 

i 



(b) Embedded crack In a plate 


T 

t 

1 


T 

2 1 

i 



(d) Corner cracks at a 
circular hole in a plate 


T 

t 

i 



T 


-- 2t 


i 


(f) Embedded cracks at a 
circular hole in a plate 


Figure 1: Crack Configurations. 

( Elliptic crack: Semi-major axis= c ; Semi-minor axis= a ) 
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Figure 2: Surface crack in a finite plate. 


2 





0>) Bending about z-axis (c) Bending about x-axis 


Oy = 1 - 2x/W 

o y =1 + 2 ( z / 1) 

0 «x < W 

0« z«-t 


Figure 3: Remote loading applied to the models. 
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44 


(NLAYER-8, NSINQU«8) 





b«~l 
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I | k I m n p q 


Singularity 

Elamant 


I J k I I J m n 



Definitions for the six faces of the Hex 8 element 


IFACE 



l-J-m-n £«0) 
l-q-p-k £=1 ) 
k-p-n-J (t|>0 ) 
l-m-q-l (q«1 ) 
H-k-J (C=0) 
p-q-m-n (C=l) 



Figure. S: Definition of Hex-8 and singularity elements. 
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(a) Twitted Element. 



% 


(b) Improperly defined element. 

( I J k I m n p q ) 

Figure. 6 : (consistently defined elements. 

( Zero or negative volumes will result for these elements.) 




(T) Element Number 
i Node Number 



KELEM(i f 1)s J 
KELEM(i,2)= K 
KELEM(i,3)= L 
KELEM(I,4)= M 
KELEM(I,5)= N 


Figure 7: Definitions and various Input parameters for the force method 







Crack front 



Forces Used In the Force Method 


FCENT(i,1,1,IR)= F ki 

FTIP(l,1,IR) = lj 

FCENT(i,1,2,IR) = ^ 

FCENT(l f 2,1,IR) = E 
■l 

FCENT(I,2,2,IR) = F 

FTIP(l,2,IR) = F b 

A 

« ft 

1= Layer number 

m ft 

IR= Current right 

ft ft 

FCENT(I,10,1,IR)= F 
h l 

FCENT(I,10,2,IR) = F h 

hand side 


COD Nodes 

IC0D(i,1) * k' 
IC0D(i,2) = I* 

IC0D(I,3) = m' 
IC0D(i,4) = n 1 
IC0D(I,5) = p* 


Figure 8: Forces and displacements in the 


layer 
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f MAXKE-5, NSINGU=8 

A 

(KELEM(I.I), U1.MAXKE) * 1 9 17 25 

33 

(MNOOE(1 , 1), 1=1 , 2*MAXKE) = 2 11 20 29 

38 

153 162 171 180 

189 

(MTIP(1,I), 1=1,2)= 1 152 

(NTIP(1,I),U1,NSINGU)= 1 2 3 4 5 

6 7 8 

(ICOD(1,l), 1=1, MAXKE) = 10 19 28 37 46 


J 


Figure 9: Example of a base model with 151 nodes and 128 elements 
(details of data input for the first layer and only 
details near the crack front are shown) 
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N LAYERs 4 NSIFsNLAYER+1 =5 
MAXKEs 5 NSINGU=8 


KELEM(NSIF,MAXKE) = 1 9 17 25 33 

129 137 145 153 161 
357 365 373 381 389 
485 493 501 509 517 

MNODE( NSIF, 2*MAXKE) = 2 11 20 29 38 

153 162 171 180 189 

153 162 171 180 189 
304 313 322 331 340 

304 313 322 331 340 
455 464 473 482 491 

455 464 473 482 491 
606 615 624 633 642 

MTIP ( NSIF, 2) s 1 152 
152 303 
303 454 
454 605 


NTI P(NSIF.NSINGU) s 1 2 3 4 5 6 7 8 

129 130 131 132 133 134 135 136 
357 358 359 360 361 362 363 364 
485 486 487 488 489 490 491 492 


ICOD(NLAYER+1 ,M AXKE) & 10 19 28 37 46 

161 170 179 188 197 

312 321 330 339 348 

463 472 481 490 499 

614 623 632 641 650 

Figure 10: (Concluded.) 
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Table 1: Input file dexl for Example 1 


SURFACE CRACK IN A PLATE TENSION AND BENDING 
SHORT 

0.30000E+08 0.30000E+00 

2161 1664 

1 

1.000000000 

0.000000000 

2 

1.013200000 

0.000000000 

3 

1.012200000 

0.005200000 

4 

1.009300000 

0.009300000 

5 

1.005200000 

0.012200000 

6 

1.000000000 

0.013200000 

7 

0.994800000 

0.012200000 

8 

0.990700000 

0.009300000 

9 

0.987800000 

0.005200000 

10 

0.986800000 

0.000000000 


A/C 


2161 

25. 

,000000000 


125.000000000 


1 210 

1 

2 

211 

210 

1 3 

212 

1 

2 210 

1 

3 

212 

210 

1 4 

213 

1 

3 210 

1 

4 

213 

210 

1 5 

214 

1 

4 210 

1 

5 

214 

210 

1 6 

215 

1 

1663 2146 

2145 

2159 

2160 

2090 

2089 2103 

2104 

0 

1664 2147 

2146 

2160 

2161 

2091 

2090 2104 

2105 

0 

1 

0 

1 

0 





; 2 

0 

1 

0 





11 
« 1 • 

0 

1 

0 





* « • 
1878* 

1 

0 

0 





1879 

1 

0 

0 





1880 

1 

0 

0 





1881 

1 

0 

0 





2161 

0 

0 

1 





0 

0 

0 

0 





3 








REMOTE 








177 

177 

1 

1 

4 




178 

178 

1 

1 

4 




179 

179 

1 

1 

4 




1625 

1625 

1 

1 

4 




1638 

1638 

1 

1 

4 




1651 

1651 

1 

1 

4 




1664 

1664 

1 

1 

4 




0 

0 

0 

0 

0 





1.0 A/T-0.2 


0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 

0.000000000 


-5.000000000 
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203 

0.0000 

1.0000 

0.0000 

204 

0.0000 

1.0000 

0.0000 

205 

0.0000 

1.0000 

0.0000 

206 

0.0000 

1.0000 

0.0000 

207 

0.0000 

1.0000 

0.0000 


2147 

0.0000 

1.0000 

0.0000 

2161 

0.0000 

1.0000 

0.0000 

0 

0.0000 

0.0000 

0.0000 

203 

0.0000 

1.0000 

0.0000 

204 

0.0000 

1.0000 

0.0000 

205 
• • » 

0.0000 

1.0000 

0.0000 

• • * 
2i05 

0.0000 

-0.5000 

0.0000 

2119 

0.0000 

- 1.0000 

0.0000 

2133 

0.0000 

- 1.0000 

0.0000 

2147 

0.0000 

- 1.0000 

0.0000 

2161 

0.0000 

- 1.0000 

0.0000 

0 

0.0000 

0.0000 

0.0000 

203 

0.0000 

0.6000 

0.0000 

204 

0.0000 

0.6800 

0.0000 


2091 

0.0000 

-0.3600 

0.0000 


2105 

0.0000 

- 1.0000 

0.0000 


2119 

0.0000 

0.4400 

0.0000 


2133 

0.0000 

0.1200 

0.0000 


2147 

0.0000 

-0.3600 

0.0000 


2161 

0.0000 

- 1.0000 

0.0000 


0 

0.0000 

0.0000 

0.0000 


0 

1 

0 0 

0 0 . 00000E+00 

0. OOOOOE+OO 

1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


1310 

1311 

1312 

1313 

1314 

1315 

1316 

1317 

1318 

1319 

1320 

1321 

1322 

1323 

1324 

1325 

1 





1 

0.0000 





8 8 





10 

19 

28 

37 

46 

219 

228 

237 

246 

255 

428 

437 

446 

455 

464 

637 

646 

655 

664 

673 


0 . OOOOOE+OO 
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846 

855 

1055 

1064 

1264 

1273 

1473 

1482 

1682 

1691 

1 

210 

628 

628 

1046 

1255 

1673 


1 

2 

6 

7 

185 

186 

190 

365 

369 

370 

548 

549 

553 

554 

732 

733 

911 

912 

916 

917 

1095 

1096 

1100 

1275 

1279 

1280 

1 

9 

183 

191 

365 

373 

547 

555 

729 

737 

911 

919 

1093 

1101 

1275 

1283 

2 

11 

211 

220 

211 

220 

420 

429 

420 

429 

629 

638 

629 

638 

838 

847 

838 

847 

1047 

1056 

1047 

1056 

1256 

1265 

1256 

1265 

1465 

1474 

1465 

1474 

1674 

1683 

125.0000 

25.0000 


864 

873 

1073 

1082 

1282 

1291 

1491 

1500 

1700 

1709 

210 

419 

837 

837 

1255 

1464 

3 

4 

8 

183 

187 

188 

366 

367 

371 

372 

550 

551 

729 

730 

734 

735 

913 

914 

918 

1093 

1097 

1098 

1276 

1277 

1281 

1282 

17 

25 

199 

207 

381 

389 

563 

571 

745 

753 

927 

935 

1109 

1117 

1291 

1299 

20 

29 

229 

238 

229 

238 

438 

447 

438 

447 

647 

656 

647 

656 

856 

865 

856 

865 

1065 

1074 

1065 

1074 

1274 

1283 

1274 

1283 

1483 

1492 

1483 

1492 

1692 

1701 

0.2000 

5.0000 


882 

1091 

1300 

1509 

1718 

419 

1046 

1464 

5 

184 

189 

368 

547 

552 

731 

736 

915 

1094 

1099 

1278 

33 

215 

397 

579 

761 

943 

1125 

1307 

38 

247 

247 

456 

456 

665 

665 

874 

874 

1083 

1083 

1292 

1292 

1501 

1501 

1710 

1.0000 
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Table 2: Output file out 12 for Example 1. 


** ************************************************ # ^^^ 
SURFACE CRACK IN A PLATE TENSION AND BENDING A/C-l 0 A/T-0 ? 

♦a*********************************************^^^/^^* 


DESCRIPTION OF THE MODEL 


OUTPUT OPTION 
YOUNG S MODULUS 
POISSION S RATIO 

NUMBER OF NODES IN THE MODEL - 
NUMBER OF ELEMENTS IN THE MODEL- 

SHORT 

0.300000E+08 

0.300 

2161 

1664 

NODE 

NODAL COORDINATES 

X-COORD Y-COORD Z-COORD 


1 1.00000 

2 1.01320 

3 1.01220 

4 1.00930 


0.00000 

0.00000 

0.00520 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN- 0 


♦i*********************************************^^^^^^^^ 

LOADING NUMBER 1 

*************************************************^^^^^^^ 

******************************************* * ************************** 
EQUILIBRIUM CHECKS 
SUM OF THE X FORCE* 0. 1218192E-09 

SUM OF THE Y FORCE* -0. 25381 15E-08 

SUM OF THE Z FORCE- 0. 1493786E-09 


APPLIED LOAD AND THE 
X- COMPONENTS: FORCE- 

Y-COMPONENTS: FORCE* 

Z-COMPONENTS: FORCE- 


SURFACE AREA COMPONENTS 
0 . OOOOOOE+OO AREA- 

0.1 25000E+03 AREA* 

0 . OOOOOOE+OO AREA- 


0. OOOOOOE+OO 
0. 125000E+03 
0. OOOOOOE+OO 


************************************************^^^^^^^^ 


NOMINAL STRESSES 
NOMINAL STRESS IN THE X- DIRECTION * 

NOMINAL STRESS IN THE Y- DIRECTION * 

NOMINAL STRESS IN THE Z- DIRECTION = 


O.OOOOOOOE+OO 
0. 1000000E+01 
O.OOOOOOOE+OO 
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************************************************************ 
STRESS INTENSITY FACTORS ARE AS FOLLOWS 

********************************************************* 


STATION 

PHI 

FROM THE FORCE METHOD 
ABSOLUTE-K 

K/(S*SQRT(PI*A/Q) ) 

1 

0.000 

0.1297875E+01 

0.1149418E+01 

2 

11.250 

0. 1261908E+01 

0.11 17566E+01 

3 

22.500 

0.1215839E+01 

0. 1076766E+01 

4 

33.750 

0. 1188309E+01 

0.1052385E+01 

5 

45.000 

0 . 1171 463 E+0 1 

0.1037467E+01 

6 

56.250 

0. 1161649E+01 

0.1028774E+01 

7 

67.500 

0. 1156257E+01 

0.1024000E+01 

8 

78.750 

0.11 53536E+01 

0 . 1021 590E+01 

9 

90.000 

0.11 52694E+01 

0.1020844E+01 

1 

FROM THE 
0.000 

CRACK OPENING DISPLACEMENT METHOD 

0. 1335318E+01 0.1182579E+01 

2 

11.250 

0. 1243665E+01 

0.1101410E+01 

3 

22.500 

0.1202907E+01 

0.1065313E+01 

4 

33.750 

0.1176129E+01 

0. 1041599E+01 

5 

45.000 

0.1159989E+01 

0.1027305E+01 

6 

56.250 

0. 1150607E+01 

0.1018996E+01 

7 

67.500 

0.1145500E+01 

0. 1014473E+01 

8 

78.750 

0. 1142883E+01 

0 . 10121 55E+01 

9 

90.000 

0. 1142066E+01 

0.1011432E+01 


******************* X1 *~~~" " 
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******************************** lfcl(r ^^ #A ^^ #A# ^^^^^^^^ 

K/( S*SQRT(PI A/Q) ) 
FORCE-METHOD COD METHOD 

**********************i** 


STATION PHI 

a****************************^^”^ 


1 

2 

3 

4 

5 

6 

7 

8 
9 


0.000 

0.11494E+01 

11.250 

0.11176E+01 

22.500 

0.10768E+01 

33.750 

0.10524E+01 

45.000 

0.10375E+01 

56.250 

0.10288E+01 

67.500 

0.10240E+01 

78.750 

0.10216E+01 

90.000 

0. 10208E+01 


0. 11826E+01 
0. 11014E+01 
0. 10653E+01 
0. 10416E+01 
0. 10273E+01 
0. 10190E+01 
0.10145E+01 
0. 10122E+01 
0. 10114E+01 


*******************************^^^^^^^^^ 

_ LOADING NUMBER 2 

********************************^^ # ^^^ # ^^^^^ 


****************i t *-t,i l i t1tit1tiri , ieieie1eificil 

EQUILIBRIUM 
SUM OF THE X FORCE- 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


************** i *********************** 

CHECKS 
0.8234746E-08 
-0.4588969E-09 
-0. 1175395E-06 


X-COMPONENTS 
Y- COMPONENTS 
Z- COMPONENTS 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 


FORCE- 

FORCE 

FORCE 


0 . OOOOOOE+OO 
0.939978E-04 
0. OOOOOOE+OO 


AREA- 

AREA- 

AREA- 


0. OOOOOOE+OO 
0.125000E+03 
0. OOOOOOE+OO 


nominal stresses 

NOMINAL STRESS IN THE X- DIRECTION - 0 OOOOOOnFj.nn 

NOMINAL STRESS IN THE Y- DIRECT ON - OtSStFS 

*»*«*- STRESS IN THE Z- DIREcilSi - O.OOOOOOOMO 


******************************^^^^^^^^^^ 

STRESS INTENSITY FACTORS ARE AS FOLLOWS 
~**«**^****^^ 
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FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE -K 

K/ (S*SQRT ( PI*A/Q) ) 

1 

0.000 

0.1160053E+01 

0. 1027362E+01 

2 

11.250 

0. 1097790E+01 

0.9722205E+00 

3 

22.500 

0. 1012613E+01 

0.8967860E+00 

4 

33.750 

0.9453700E+00 

0.8372347E+00 

5 

45.000 

0.8916253E+00 

0.7896376E+00 

6 

56.250 

0.8502292E+00 

0.7529765E+00 

7 

67.500 

0. 8206783 E+00 

0. 7268058E+00 

8 

78.750 

0. 80284 58E+00 

0.7110130E+00 

9 

90.000 

0.79687 13 E+00 

0.7057219E+00 


FROM THE 

CRACK OPENING DISPLACEMENT METHOD 

1 

0.000 

0. 1208218E+01 

0. 1070017E+01 

2 

11.250 

0. 1081892E+01 

0.9581404E+00 

3 

22.500 

0. 1001718E+01 

0.8871372E+00 

4 

33.750 

0.9348640E+00 

0.8279304E+00 

5 

45.000 

0.8816859E+00 

0.7808351 E+00 

6 

56.250 

0.8405767E+00 

0.7444282E+00 

7 

67.500 

0.81 12083 E+00 

0.7184190E+00 

8 

78.750 

0.7934419E+00 

0.7026848E+00 

9 

90.000 

0.7874846E+00 

0.6974089E+00 


************************************************************ 

NUMBER OF SINGULARITY ELEMENTS IN THE MODEL- 64 
************************************************************ 

************************************************************ 




K/( S*SQRT(PI A/Q) 

) 

STATION 

PHI 

FORCE -METHOD 

COD METHOD 

************************************************************ 

1 

0.000 

0.10274E+01 

0. 10700E+01 

2 

11.250 

0.97222E+00 

0.95814E+00 
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3 

22.500 

0.89679E+00 

0.88714E+00 

4 

33.750 

0.83723E+00 

0.82793E+00 

5 

45.000 

0.78964E+00 

0.78084E+00 

6 

56.250 

0.75298E+00 

0.74443E+00 

7 

67.500 

0. 72681 E+00 

0.71842E+00 

8 

78.750 

0.71101E+00 

0.70268E+00 

9 

90.000 

0.70572E+00 

0.69741 E+00 


********************************************************************** 

LOADING NUMBER 3 

********************************************************************** 


********************************************************************** 


EQUILIBRIUM 
SUM OF THE X FORCE- 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


CHECKS 
0.1084288E-10 
-0.3993250E-10 
-0 . 471 1880E- 10 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 
X-COMPONENTS: FORCE- O.OOOOOOE+OO AREA- O.OOOOOOE+OO 

Y-COMPONENTS: FORCE- 0.251484E-04 AREA- 0.125000E+03 

Z-COMPONENTS: FORCE- O.OOOOOOE+OO AREA- O.OOOOOOE+OO 


A********************************************************************* 

NOMINAL STRESSES 
NOMINAL STRESS IN THE X- DIRECTION - 0 . OOOOOOOE+OO 

NOMINAL STRESS IN THE Y- DIRECTION - 0 . 201 1872E-06 

NOMINAL STRESS IN THE Z- DIRECTION - 0. OOOOOOOE+OO 


************************************************************ 


STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************************************************************ 


STATION 

PHI 

1 

0.000 

2 

11.250 

3 

22.500 

4 

33.750 


FROM THE FORCE METHOD 
ABSOLUTE-K 
0.2430848E-02 
0.2362947E-02 
0.2276044E-02 
0.2224167E-02 


K/(S*SQRT(PI*A/Q) ) 
0.2152797E-02 
0.2092664E-02 
0.2015700E-02 
0. 1969758E-02 
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5 

45.000 

0.2192557E-02 

0. 1941763E-02 

6 

56.250 

0.2174276E-02 

0.1925574E-02 

7 

67.500 

0.2164344E-02 

0. 1916778E-02 

8 

78.750 

0 . 2159391 E-02 

0 . 1912391 E-02 

9 

90.000 

0.2157872E-02 

0.1911046E-02 


FROM THE 

CRACK OPENING DISPLACEMENT METHOD 

1 

0.000 

0.2501464E-02 

0.2215336E-02 

2 

11.250 

0.2328975E-02 

0.2062578E-02 

3 

22.500 

0.2251980E-02 

0. 1994389E-02 

4 

33.750 

0.2201486E-02 

0 . 1949671 E-02 

5 

45.000 

0.2171 161E-02 

0. 1922815E-02 

6 

56.250 

0.2153651E-02 

0. 1907308E-02 

7 

67.500 

0.2144214E-02 

0. 1898950E-02 

8 

78.750 

0.2139433E-02 

0. 1894716E-02 

9 

90.000 

0.2137952E-02 

0. 1893404E-02 


*********** r************************************************* 


NUMBER OF SINGULARITY ELEMENTS IN THE MODEL* 64 
************************************************************ 


************************************************************ 


K/( S*SQRT(PI A/Q) ) 

STATION PHI FORCE-METHOD COD METHOD 

************************************************************ 


1 

0.000 0.21528E-02 

0.22153E-02 


2 

11.250 0.20927E-02 

0.20626E-02 


3 

22.500 0.20157E-02 

0. 19944E-02 


4 

33.750 0. 19698E-02 

0. 19497E-02 


5 

45.000 0. 19418E-02 

0. 19228E-02 


6 

56.250 0. 19256E-02 

0. 19073E-02 


7 

67.500 0.19168E-02 

0.18989E-02 


8 

78.750 0.19124E-02 

0. 18947E-02 
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9 


90.000 


0.19110E-02 


0.18934E-02 


ALL ELEMENTS SATISFY EQUILIBRIUM 
STOP (called by $MAIN ) 

CP: 29.608s, Wall clock: 84.966s, 8.7% of 4-CPU Machine 
HUM mem: 7733236, HWM stack: 310499, Stack overflows: 0 


Table 3: Output file outn!2 for Example 2. 


**********************************************************^ 
SURFACE CRACK IN A PLATE TENSION AND BENDING A/C-1.0 A/T-0 2 
**********************************************************;^ w * 


DESCRIPTION OF THE MODEL 


OUTPUT OPTION 
YOUNG S MODULUS 
POISSION S RATIO 

NUMBER OF NODES IN THE MODEL - 
NUMBER OF ELEMENTS IN THE MODEL- 


NODAL COORDINATES 

NODE X-COORD Y-COORD Z-COORD 


SHORT 

0.300000E+08 

0.300 

2161 

1664 


1.00000 0.00000 0.00000 

1.01320 0.00000 0.00000 

1-01220 0.00520 0.00000 


***********************************1c**1'1'1,1'1c1 [ *- k 1'1'1'1'1'i'1'1' i '1' i ' i '1' 1 ' ic i' 1 ' i ' iH ' i ' 1 ' il ' 1 ' 


equilibrium 

SUM OF THE X FORCE- 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


CHECKS 
0.1052262E-09 
-0.2110426E-08 
0. 1165562E-09 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 
X-COMPONENTS: FORCE- 0 . 000000E+00 AREA- 

Y-COMPONENTS: FORCE- 0.125000E+03 AREA- 

Z-COMPONENTS: FORCE- 0.000000E+00 AREA- 


0 . OOOOOOE+OO 
0.125000E+03 
0. OOOOOOE+OO 
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^************************************************** 
****************** | ctdf^^FS 

N c?c M tn THF X DIRECTION - 0.0000000E+00 

nominal stress in THE X- DIRtCT ON _ 0.1000000E+01 

Z\Z STRESS » m i- DIRECTION O.OOOOOOOE+OO 

***************************************************** 

* 5 ********************************************************** 

^^******************************************** 

*********** K/( S*SQRT(PI VQ) ) 

nni FROM VCCT FROM VCCT 

STATION p-STRAIN P-STRESS 

r J I ivun . 


1 

2 

3 

4 

5 

6 

7 

8 
Q 


0.000 

0. 11774E+01 

0.11232E+01 

11.250 

0. 11215E+01 

0 . 10699E+01 

22.500 

0. 10819E+01 

0 . 10320E+01 

33.750 

0.10576E+01 

0. 10089E+01 

45.000 

0.10428E+01 

0.99475E+00 

56.250 

0.10341E+01 

0.98650E+00 

67.500 

0.10294E+01 

0.98196E+00 

78.750 

0.10269E+01 

0.97964E+00 

90.000 

0.10262E+01 

0.97892E+00 


S2K2SJ5S5S5 — — 

rrrrrrrrrrrrr^ 

SS OF T THE ? force'- 

SUM OF THE Z FORCE- -0.985Z94UL u/ 


applied LOAD AND 
X-COMPONENTS: FORCE= 
Y- COMPONENTS: 

Z- COMPONENTS: 


FORCE" 

FORCE* 


THE SURFACE AREA COMPONENTS 
0.000000E+00 AREA- 

0.939978E-04 AREA- 

0 000000E+00 area= 


O.OOOOOOE+OO 
0.125000E+03 
0 . 000000E+00 
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***************itlt*le1t*1t 


*********+***l,iH'+ i m 1til ', 

nominal STOP 
SK 25™ » !(!§ x- DIRECTION 

direction - 

DIRECTION - 


^ aii int a 

NOMINAL STRESS IN THE Y 
NOMINAL STRESS IN THE Z 


************************** 
S S E S 

0.0000000E+00 
0.7519827E-06 
O.OOOOOOOE+OO 


*********************************^JLjt 4 ^ . . 

NUMBER OF SINGULARITY ELEMFNTC tm tuc ******************** 

................ 

STRTION PHI F &$ «}"<" VQ) ) 


1 

2 

3 

4 

5 

6 

7 

8 
9 


0.000 

0.10588E+01 

11.250 

0.97573E+00 

22.500 

0.90118E+00 

33.750 

0.84131E+00 

45.000 

0.79357E+00 

56.250 

0.75670E+00 

67.500 

0.73037E+00 

78.750 

0.71446E+00 

90.000 

0.70912E+00 


0.10100E+01 

0.93078E+00 

0.85967E+00 

0.80256E+00 

0.75701E+00 

0. 72185E+00 

0.69673E+00 

0.68155E+00 

0.67646E+00 


..................... — 

e q u i l i*B*R*ru*M***c*H*rrrr******”*********“*'** 

SUM OF THE Y FORCE* 0.1055867E-10 

aSSJaS :S:K:JX 


X-COMPONENTS 

Y-COMPONENTS 

Z-COMPONENTS 


APPLIED LOAD AND 


FORCE- 

FORCE- 

FORCE- 


THE SURFACE AREA 
O.OOOOOOE+OO 
0.251484E-04 
O.OOOOOOE+OO 


COMPONENTS 

AREA* 

AREA* 

AREA* 


O.OOOOOOE+OO 

0.125000E+03 

O.OOOOOOE+OO 


64 



****************************************** ## ^^ # ^ lk ^ A ^^^^^ 

NOMINAL STRESSES 
NOMINAL STRESS IN THE X- DIRECTION - 0 OOOOOOOF + nn 

NOMINAL STRESS IN THE Y- DIRECTION - 0 2011872F M 

NOMINAL STRESS IN THE Z- DIRECTION - O.'OOOOOOOE+OO 


***********************************^^^^^^^ 
NUMBER OF SINGULARITY ELEMENTS IN THE MODEL* 0 

********************************** 1 , 1 ,^^^^^^^^ 

********************************** #a ^ aa#a#aa# ^ #### ^ #### ^^^ 
STATION PH, FROM Tc^ W > l m VCCT 

*******************************J;^^^ a ^^ a ^p-stress 


1 

0.000 

0.22053E-02 

0.21038E-02 

2 

11.250 

0.21001E-02 

0.20033E-02 

3 

22.500 

0.20253E-02 

0. 19320E-02 

4 

33.750 

0. 19795E-02 

0. 18883E-02 

5 

45.000 

0. 19517E-02 

0. 18618E-02 

6 

56.250 

0.19356E-02 

0. 18465E-02 

7 

67.500 

0.19269E-02 

0.18382E-02 

8 

78.750 

0. 19225E-02 

0. 18339E-02 

9 

90.000 

0. 19211E-02 

0. 18327E-02 


ALL ELEMENTS SATISFY EQUILIBRIUM 


STOP (called by SMAIN ) 

CP: 30.631s, Wallclock: 37.531s, 20.4% of 4-CPU Machine 

HWM mem: 7730259, HWM stack: 310502, Stack overflows: 0 


Table 4: Input file dex3a for Example 3. 


SURFACE CRACK IN A PLATE 
SHORT 

0.30000E+08 0.30000E+00 

2441 1872 

1 5.000000000 

2 5.002656718 

3 5.002454282 

4 5.001868300 


A/C=0.2 A/T=0.2 


0.000000000 

0.000000000 

0.005200000 

0.009300000 


0.000000000 

0.000000000 

0.000000000 

0.000000000 
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» • • 


2437 25.000000000 

2438 25.000000000 

2439 25.000000000 

2440 25.000000000 

2441 25.000000000 

1 210 1 2 211 

2 210 1 3 212 

3 210 1 4 213 

4 210 1 5 214 

5 210 1 6 215 

6 210 1 7 216 



15.000000000 



25.000000000 



45.000000000 



85.000000000 



125.000000000 


210 

1 

3 

212 

1 

210 

1 

4 

213 

1 

210 

1 

5 

214 

1 

210 

1 

6 

215 

1 

210 

1 

7 

216 

1 

210 

1 

8 

217 

1 


1867 2422 

1868 2423 

1869 2424 

1870 2425 

1871 2426 

1872 2427 

1 

2 

11 


2421 2435 

2422 2436 

2423 2437 

2424 2438 

2425 2439 

2426 2440 
0 1 


2436 2310 

2437 2311 

2438 2312 

2439 2313 

2440 2314 

2441 2315 
0 


2309 2323 2324 

2310 2324 2325 

2311 2325 2326 

2312 2326 2327 

2313 2327 2328 

2314 2328 2329 


1 

1 


0 

0 

0 

0 

0 

0 


1878 1 0 0 

1879 1 0 0 

1880 1 0 0 

1881 1 0 0 

2441 0 0 1 

0 0 0 0 

1 

REMOTE 

177 177 1 1 

178 178 1 1 

179 179 1 1 


4 

4 

4 


1859 1859 1 

1872 1872 1 

0 0 0 

203 0.0000 

204 0.0000 

205 0.0000 
206 0.0000 


2399 0.0000 
2413 0.0000 
2427 0.0000 
2441 0.0000 


1 4 

1 4 

0 0 

1.0000 0.0000 

1.0000 0.0000 

1.0000 0.0000 

1.0000 0.0000 


1.0000 0.0000 

1.0000 0.0000 

1.0000 0.0000 

1.0000 0.0000 


- 5.000000000 

- 5.000000000 

- 5.000000000 

- 5.000000000 

- 5.000000000 
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0 . OOOOOE+OO O.OOOOOE+OO 


0 

0 

1 

1 

6 

11 

16 


0.0000 

0 0 


2 

7 

12 

17 


0.0000 0.0000 

0 0.00000E+00 

3 4 
8 9 
13 14 
18 19 


1590 

1591 

1595 

1596 

1600 

1601 

1605 


1 


1 


.0000 


8 


10 

19 

219 

228 

428 

437 

637 

646 

846 

855 

1055 

1064 

1264 

1273 

1473 

1482 

1682 

1691 

1 

210 

628 

628 

1046 

1255 

1673 


1 

2 

6 

7 

185 

186 

190 

365 

369 

370 

548 

549 

553 

554 

732 

733 

911 

912 

916 

917 

1095 

1096 

1100 

1275 

1279 

1280 

1 

9 

183 

191 

365 

373 

547 

555 

729 

737 

911 

919 

1093 

1101 

1275 

1283 

2 

11 

211 

220 

211 

220 

420 

429 


1592 1593 
1597 1598 
1602 1603 


28 

37 

237 

246 

446 

455 

655 

664 

864 

873 

1073 

1082 

1282 

1291 

1491 

1500 

1700 

1709 

210 

419 

837 

837 

1255 

1464 

3 

4 

8 

183 

187 

188 

366 

367 

371 

372 

550 

551 

729 

730 

734 

735 

913 

914 

918 

1093 

1097 

1098 

1276 

1277 

1281 

1282 

17 

25 

199 

207 

381 

389 

563 

571 

745 

753 

927 

935 

1109 

1117 

1291 

1299 

20 

29 

229 

238 

229 

238 

438 

447 


5 

10 

15 

20 


1594 

1599 

1604 


46 

255 

464 

673 

882 

1091 

1300 

1509 

1718 

419 

1046 

1464 

5 

184 

189 

368 

547 

552 

731 

736 

915 

1094 

1099 

1278 

33 

215 

397 

579 

761 

943 

1125 

1307 

38 

247 

247 

456 
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420 

429 

629 

638 

629 

638 

838 

847 

838 

847 

1047 

1056 

1047 

1056 

1256 

1265 

1256 

1265 

1465 

1474 

1465 

1474 

1674 

1683 

125.0000 

25.0000 


438 

647 

647 

856 

856 

1065 

1065 

1274 

1274 

1483 

1483 

1692 

0.2000 


447 

656 

656 

865 

865 

1074 

1074 

1283 

1283 

1492 

1492 

1701 

5.0000 


456 

665 

665 

874 

874 

1083 

1083 

1292 

1292 

1501 

1501 

1710 

0.2000 


Table 5: Output file outr22 for Example 3(a). 


******************************** <r#ifc## ^^ ##A# ^ <k ^^^ 

SURFACE CRACK IN A PLATE A/C-0.2 A/T-0 2 
*******************************^^^^VT # 0.2^ 



DESCRIPTION OF THE MODEL 



OUTPUT OPTION 
YOUNG S MODULUS 
POISSION S RATIO 

NUMBER OF NODES IN THE MODEL - 
NUMBER OF ELEMENTS IN THE MODEL- 

SHORT 

0.300000E+08 

0.300 

2441 

1872 

NODE 

NODAL COORDINATES 

X- COORD Y-COORD Z-COORD 


1 

2 

3 

4 


5.00000 

5.00266 

5.00245 

5.00187 


0.00000 

0.00000 

0.00520 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN* 


LOADING NUMBER 1 

*****************************^^^^^^^^^^^^^^^^^^^ 

EQUILIBRIUM CHECKS 
SUM OF THE X FORCE* 0.4193217E-09 

SUM OF THE Y FORCE* -0.3798505E-08 

SUM OF THE Z FORCE- -0.4241372E-09 
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X- COMPONENTS 
Y-COMPONENTS 
Z- COMPONENTS 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 


FORCE* O.OOOOOOE+OO AREA- 
FORCE- 0. 124996E+03 AREA- 
FORCE- O.OOOOOOE+OO AREA- 


O.OOOOOOE+OO 

0.124996E+03 

O.OOOOOOE+OO 


NOMINAL STRESSES 
NOMINAL STRESS IN THE X- DIRECTION - 0 OOOOOOOF^nn 

NOMINAL STRESS IN THE Y- DIRECTION - o’iOOOoSeU? 

N « INA L stress in the z- direction - o:iooooooMo 


**“**™*******™******™******1M***1«M 1Ht1t1ti ' 1 ' 1l1t1 ' 1 ' i ' 1 „ 1 ' i ' 1k1 ' 

. . o STRESS INTENSITY FACTORS ARE AS FOLLOWS 
******************************** a ^^^^lluw^^^^^^^ 


FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE-K 

K/(S*SQRT(PI*A/Q)) 

1 

0.000 

0. 1042638E+01 

0.6177579E+00 

2 

11.250 

0. 1097014E+01 

0.6499757E+00 

3 

22.500 

0. 1265662E+01 

0.7498990E+00 

4 

33.750 

0 . 1461 145E+01 

0.8657214E+00 

5 

45.000 

0. 1629536E+01 

0.9654924E+00 

6 

56.250 

0. 1761401E+01 

0. 1043622E+01 

7 

67.500 

0. 1856268E+01 

0.1 09983 1E+01 

8 

78.750 

0.19 1 3422E+01 

0. 1133694E+01 

9 

90.000 

0. 1932516E+01 

0. 1145007E+01 


FROM THE 

CRACK OPENING DISPLACEMENT 

METHOD 

1 

0.000 

0.8606799E+00 

0.5099489E+00 

2 

11.250 

0. 1028380E+01 

0.6093103E+00 

3 

22.500 

0. 1265637E+01 

0. 7498842E+00 

4 

33.750 

0. 1475477E+01 

0.8742133E+00 

5 

45.000 

0. 1642917E+01 

0.9734206E+00 

6 

56.250 

0. 1769638E+01 

0. 1048502E+01 

7 

67.500 

0.1859042E+01 

0. 1101474E+01 
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8 

9 


78.750 

90.000 


0. 1912002E+01 
0.1929571E+01 


0.1132852E+01 

0.1143262E+01 


************************************************************ 
NUMBER OF SINGULARITY ELEMENTS IN THE MODEL* 64 
************************************************************ 


************************************************************ 

K/( S*SQRT(PI A/Q) ) 

station PHI FORCE-METHOD COD METHOD 

************************************************************ 

1 

0.000 

0.61776E+00 

0.50995E+00 

2 

11.250 

0.64998E+00 

0.60931E+00 

3 

22.500 

0.74990E+00 

0.74988E+00 

4 

33.750 

0.86572E+00 

0.87421E+00 

5 

45.000 

0 . 96549E+00 

0.97342E+00 

6 

56.250 

0. 10436E+01 

0. 10485E+01 

7 

67.500 

0.10998E+01 

0. 11015E+01 

8 

78.750 

0.11337E+01 

0.1 1329E+01 

9 

90.000 

0. 11450E+01 

0.11433E+01 


ALL ELEMENTS SATISFY EQUILIBRIUM 

^ P 38 ( 7&: d SlStaki 92.835s, .0.4* of 4-CPU Machine 

HUM mem: 7733236, HWM stack: 310499, Stack overflows: 0 


Table 6: Output file outc22 for Example 3(b). 


************************************************************ 
SURFACE CRACK -CRACK FACE PRESSURE LOADING A/C-0.2 A/T-0.2 
************************************************************ 


DESCRIPTION OF THE MODEL 


OUTPUT OPTION 

YOUNG S MODULUS 

POISSION S RATIO 

NUMBER OF NODES IN THE MODEL 

NUMBER OF ELEMENTS IN THE MODEL 


SHORT 

0 . 300000E+08 
0.300 


2441 

1872 
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NODE 


NODAL COORDINATES 
X-COORD Y-COORD 


Z- COORD 


1 

2 

3 

4 


5.00000 

5.00266 

5.00245 

5.00187 


0.00000 

0.00000 

0.00520 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN- 0 

********************************************************************** 

LOADING NUMBER 1 

*************************************************************** 

********************************************************************** 
equilibrium checks 
SUM OF THE X FORCE- 0. 1027622E-11 

SUM OF THE Y FORCE- -0.9257040E-11 

SUM OF THE Z FORCE- -0.5503000E-09 


X 

Y 

Z 


APPLIED LOAD AND THE 
-COMPONENTS: FORCE- 

-COMPONENTS: FORCE- 

-COMPONENTS: FORCE- 


SURFACE AREA COMPONENTS 
O.OOOOOOE+OO AREA- 

0.389798E+01 AREA- 

O.OOOOOOE+OO AREA- 


0. 175788E-16 
0.389798E+01 
0.278156E-16 


******************************* 


*************************************** 


NOMINAL STRESS 
NOMINAL STRESS IN THE X- DIRECTION - 

NOMINAL STRESS IN THE Y- DIRECTION = 

NOMINAL STRESS IN THE Z- DIRECTION = 


E S 


O.OOOOOOOE+OO 

0.1000000E+01 

O.OOOOOOOE+OO 


************************************************************ 
STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************* 2 ********************************************** 


STATION PHI 

1 0.000 

2 11.250 

3 22.500 

4 33.750 

5 45.000 


FROM THE FORCE METHOD 
ABSOLUTE-K 
0. 1038225E+01 
0. 1095227E+01 
0.1265707E+01 
0. 1462723E+01 
0. 1632234E+01 


K/(S*SQRT(PI*A/Q) ) 
0.6151432E+00 
0.6489169E+00 
0.7499253E+00 
0.8666567E+00 
0.9670910E+00 
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6 

56.250 

0. 1764889E+01 

0.1045689E+01 

7 

67.500 

0. 1860285E+01 

0. 1102210E+01 

8 

78.750 

0. 1917742E+01 

0.1136254E+01 

9 

90.000 

0.1936936E+01 

0. 1147626E+01 


FROM THE 

CRACK OPENING DISPLACEMENT 

METHOD 

1 

0.000 

0.8465032E+00 

0.5015492E+00 

2 

11.250 

0. 1021579E+01 

0. 605281 1E+00 

3 

22.500 

0. 1263493E+01 

0.7486140E+00 

4 

33.750 

0. 1474976E+01 

0.8739167E+00 

5 

45.000 

0. 1643138E+01 

0.9735514E+00 

6 

56.250 

0. 1770183E+01 

0. 1048826E+01 

7 

67.500 

0. 1859740E+01 

0. 1101887E+01 

8 

78.750 

0. 1912768E+01 

0. 1133306E+01 

9 

90.000 

0.1930357E+01 

0. 1143727E+01 


********************************^^^^^^^^ 

NUMBER OF SINGULARITY ELEMENTS IN THE MODEL- 64 
*****************************^^^^^^^ 4 ^^^^ 




STATION dut K/ ( $*SQRT(PI A/Q) ) 

1 

0.000 

0.61514E+00 

0.50155E+00 

2 

11.250 

0.64892E+00 

0.60528E+00 

3 

22.500 

0.74993E+00 

0. 74861 E+00 

4 

33.750 

0.86666E+00 

0.87392E+00 

5 

45.000 

0.96709E+00 

0.97355E+00 

6 

56.250 

0. 10457E+01 

0. 10488E+01 

7 

67.500 

0. 11022E+01 

0. 11019E+01 

8 

78.750 

0.11363E+01 

0. 11333E+01 

9 

90.000 

0.11476E+01 

0. 11437E+01 
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all elements satisfy equilibrium 


S* ( ST Vallclock: 80.092s 12 

mA ^7733236, HUM stack: 310499, 


.1% of 4 -CPU Machine 
Stack overflows: 0 


Tabie 7: Output"fiie"outr28 for Example 4(a). 


******************** 
EM 

******************* 


*********************** 


SURFACE LRACK : REfWTE^TENSI0N^^^A/C-0^2^^A/T=0^8^^^^^^^^^^^ 


NODE 


description of the model 


OUTPUT OPTION 
YOUNG S MODULUS 

numb!r°of S nodes°in the model 

NUMBER OF ELEMENTS IN THE MODEL- ^185 


= SHORT 

* 0.300000E+08 
= 0.300 
» 2464 


NODAL COORDINATES 


X-COORD 


Y-COORD 


Z-COORD 


1 

2 

3 

4 

5 


5.00000 

5.00266 

5.00245 

5.00187 

5.00102 


0.00000 

0.00000 

0.00510 

0.00930 

0.01220 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN- 0 

\ « of'the'* 6 FORCE- 0 8482506E-09 

SUM OF IHt a cfi<:Q97i>F-08 

IS °o f f t t h h e e v z JS& 48SW 
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*-~n A t P s PUE0 W 

FORCE* 


Y- COMPONENTS 


AREA- 


area. 


Z-CCRONENTS: FoS£ SlXXgg Jg* [«Kg 

_ 8Kfc A= 0 . OOOOOOE+OO 

nominal ******************************* 

NOMINAL STRESS IN THE X- DIRECTION S S E S 

NOMINAL STRESS IN THE Y- DIRECTION " O.OOOOOOOE+OO 

NOMINAL STRESS IN THE Z- D D R R T T o° " 0.1000000E+01 

o.oooooooe+oo 

**********************^^^^ 


STATION PHI 

1 0.000 

2 11.250 

3 22.500 

4 33.750 

5 45.000 

6 56.250 

7 67.500 

8 78.750 

9 90.000 


FROM THE FORCE METHOD 

absolute-k 

0.1942565E+01 

0.1 9453 17E+01 

0.2148890E+01 

0.2397951E+01 

0.2637838E+01 

0.281 51 54E+01 

0.2934686E+01 

0.2980024E+01 

0.2990430E+01 


K/ (S*SQRT ( PI*A/Q ) ) 
0. 1150961E+01 
0.1152591E+01 
0.1273207E+01 
0. 1420775E+01 
0. 1562907E+01 
0. 1 667966E+01 
0. 1738788E+01 
0. 1765651E+01 
0.1771816E+01 


1 

2 

3 

4 

5 

6 
7 


FROM THE 
0.000 
11.250 
22.500 
33.750 
45.000 
56.250 
67.500 


CRACK OPENING DISPLACEMENT METHOD 

0. 1667430E+01 0.9879444E+00 
0. I837670E+0I 0 . 108881 ! E + 01 
0.2173179E+01 0.1287598E+01 
0.2450039E+ 01 0. 1451637E + 01 
0.2 625023E+01 0.1555314E+01 
0.284 1981 E+01 0. 1 683861 E+01 
0.2935053E+01 0. 1739006E+01 
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8 


78.750 


0.2960373E+01 


0. 1754007E+01 


9 90.000 0.2963058E+01 0. 1755598E+01 

************************************************************ 

NUMBER OF SINGULARITY ELEMENTS IN THE MODEL* 64 
************************************************************ 

************************************************************ 

K/( S*SQRT ( P I A/Q) ) 

STATION PHI FORCE -METHOD COD METHOD 

************************************************************ 


1 

0.000 

0.11510E+01 

0.98794E+00 

2 

11.250 

0.11 526E+01 

0. 10888E+01 

3 

22.500 

0. 12732E+01 

0. 12876E+01 

4 

33.750 

0. 14208E+01 

0. 14516E+01 

5 

45.000 

0. 15629E+01 

0. 15553E+01 

6 

56.250 

0. 16680E+01 

0. 16839E+01 

7 

67.500 

0. 17388E+01 

0. 17390E+01 

8 

78.750 

0. 17657E+01 

0. 17540E+01 

9 

90.000 

0. 17718E+01 

0. 17556E+01 


ALL ELEMENTS SATISFY EQUILIBRIUM 
STOP (called by SMAIN ) 

CP: 45.682s, Wallclock: 85.788s, 13.3% of 4-CPU Machine 

HUM mem: 7733236, HWM stack: 310499, Stack overflows: 0 


Table 8: Output file outc28 for Example 4(b). 


************************************************************ 

SURFACE CRACK-CRACK FACE PRESSURE LOADING A/C=0.2 ,A/T=0.8 
************************************************************ 


DESCRIPTION OF THE MODEL 


OUTPUT OPTION 

33 

SHORT 

YOUNG S MODULUS 

33 

0.300000E+08 

POISSION S RATIO 

S3 

0.300 

NUMBER OF NODES IN 

THE MODEL =* 

2464 

NUMBER OF ELEMENTS 

IN THE MODEL= 

1856 
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NODAL COORDINATES 

NODE X-COORD Y-COORD Z-COORD 


1 

2 

3 

4 


5.00000 

5.00266 

5.00245 

5.00187 


0.00000 

0.00000 

0.00510 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN- 0 

********************************************************************** 

LOADING NUMBER 1 

******r************f***f**r*************************f*******r*************** 


★Jr******************************************************************** 


EQUILIBRIUM 
SUM OF THE X FORCE- 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


CHECKS 
-0.4298795E-10 
-0. 7681 145E- 10 
-0. 1726783E-07 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 


X-COMPONENTS: 

FORCE- 

0 . OOOOOOE+OO 

AREA- 

0. 184609E-16 

Y-COMPONENTS: 

FORCE- 

0.389798E+01 

AREA- 

0.389798E+01 

Z- COMPONENTS: 

FORCE- 

0. OOOOOOE+OO 

AREA- 

0.390507E-16 


********************************************************************** 

NOMINAL STRESSES 
NOMINAL STRESS IN THE X- DIRECTION - 0 . OOOOOOOE+OO 

NOMINAL STRESS IN THE Y- DIRECTION - 0. 1000000E+01 

NOMINAL STRESS IN THE Z- DIRECTION - 0 . OOOOOOOE+OO 


***★ *r ****** * *r* * *r*r *r *r*r *r **r *r* *r* *** * * ** **** ********************** 

STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************************************************************ 


FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE-K 

K/ ( S*SQRT ( P I*A/Q ) ) 

1 

0.000 

0. 1939965E+01 

0.11 49420E+0 1 

2 

11.250 

0.1946502E+01 

0.1 153294E+01 

3 

22.500 

0.2152937E+01 

0.1275605E+01 

4 

33.750 

0.2403408E+01 

0.1424008E+01 

5 

45.000 

0.2629350E+01 

0.1557878E+01 
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6 


56.250 


0.2818348E+01 


7 

8 
9 


67.500 

78.750 

90.000 


0.2939771E+01 
0. 2986851 E+01 
0.2997954E+01 


0. 1669858E+01 
0. 1741801E+01 
0. 1769696E+01 
0. 1776274E+01 



FROM THE 

CRACK OPENING DISPLACEMENT METHOD 

1 

0.000 

0. 1654822E+01 

0.9804745E+00 

2 

11.250 

0 . 183181 2E+01 

0. 1085340E+01 

3 

22.500 

0.2171463E+01 

0. 1286582E+01 

4 

33.750 

0.2449731E+01 

0. 1451454E+01 

5 

45.000 

0.2625300E+01 

0. 1555478E+01 

6 

56.250 

0.2842561E+01 

0. 1684204E+01 

7 

67.500 

0.2935786E+01 

0. 1739440E+01 

8 

78.750 

0 . 2961 178E+01 

0. 1754485E+01 

9 

90.000 

0.2963886E+01 

0. 1756089E+01 


************************************************************ 


NUMBER OF SINGULARITY ELEMENTS IN THE MODEL= 64 

************************************************************ 


************************************************************ 


K/( S*SQRT (PI A/Q) ) 

STATION PHI FORCE-METHOD COO METHOD 

************************************************************ 


1 

0.000 0.11494E+01 

0.98047E+00 


2 

11.250 0.11533E+01 

0. 10853E+01 


3 

22.500 0. 12756E+01 

0. 12866E+01 


4 

33.750 0. 14240E+01 

0. 14515E+01 


5 

45.000 0. 15579E+01 

0. 15555E+01 


6 

56.250 0.16699E+01 

0. 16842E+01 


7 

67.500 0.17418E+01 

0.17394E+01 


8 

78.750 0. 17697E+01 

0. 17545E+01 


9 

90.000 0. 17763E+01 

0 . 17561 E+01 
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ALL ELEMENTS SATISFY EQUILIBRIUM 
STOP (called by $MAIN ) 

CP: 45.254s, Wallclock: 115.863s, 9.8% of 4-CPU Machine 

HUM mem: 7733236, HUM stack: 310499, Stack overflows: 0 


Table 9: Output file outcor28 for Example 5. 


************************************************************ 

CORNER CRACK-REMOTE TENSION A/C-0.2 ,A/T-0.8 
************************************************************ 


DESCRIPTION OF THE MODEL 


OUTPUT OPTION 

m 

SHORT 

YOUNG S MODULUS 

m 

0.300000E+08 

POISSION S RATIO 

B 

0.300 

NUMBER OF NODES IN 

THE MODEL - 

2464 

NUMBER OF ELEMENTS 

IN THE MODEL- 

1856 


NODAL COORDINATES 

NODE X-COORD Y-COORD Z-COORD 


1 

2 

3 

4 


5.00000 

5.00266 

5.00245 

5.00187 


0.00000 

0.00000 

0.00510 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN- 0 

********************************************************************** 

LOADING NUMBER 1 

********************************************************************** 


********************************************************************** 


EQUILIBRIUM 
SUM OF THE X FORCE= 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


CHECKS 
0. 1055787E-08 
-0.4741175E-08 
-0.2597931E-06 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 
X-COMPONENTS: FORCE- 0.000000E+00 AREA- O.OOOOOOE+OO 

Y-COMPONENTS: FORCE- 0.624962E+02 AREA- 0.624962E+02 

Z-COMPONENTS: FORCE- 0.000000E+00 AREA- 0.000000E+00 


78 


********************************************************************** 

NOMINAL STRESSES 
NOMINAL STRESS IN THE X- DIRECTION - O.OOOOOOOE+OO 

NOMINAL STRESS IN THE Y- DIRECTION - 0. lOOOOOOE+Ol 

NOMINAL STRESS IN THE Z- DIRECTION - O.OOOOOOOE+OO 


************************************************************ 

STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************************************************************ 


FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE-K 

K/(S*SQRT(PI*A/Q) ) 

1 

0.000 

0.2016223E+01 

0. 1194603E+01 

2 

11.250 

0.2019591E+01 

0. 1196598E+01 

3 

22.500 

0.2228682E+01 

0. 1320484E+01 

4 

33.750 

0.2490525E+01 

0. 1475624E+01 

5 

45.000 

0. 2750801 E+01 

0. 1629837E+01 

6 

56.250 

0.2961335E+01 

0. 1754578E+01 

7 

67.500 

0.3142617E+01 

0. 1861986E+01 

8 

78.750 

0.3339150E+01 

0. 1978431E+01 

9 

90.000 

0.3552545E+01 

0.2104866E+01 


FROM THE 

CRACK OPENING DISPLACEMENT METHOD 

1 

0.000 

0. 1740380E+01 

0. 1031 167E+01 

2 

11.250 

0. 1905171E+01 

0. 1128805E+01 

3 

22.500 

0.2252796E+01 

0 . 1334771 E+01 

4 

33.750 

0.2542841E+01 

0.1 506621 E+01 

5 

45.000 

0.2734719E+01 

0. 1620308E+01 

6 

56.250 

0.2983844E+01 

0. 1767914E+01 

7 

67.500 

0.31 297 1 7E+01 

0. 1854343E+01 

8 

78.750 

0.3271988E+01 

0. 1938638E+01 

9 

90.000 

0.3648129E+01 

0.2161499E+01 
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***************************************** # ** ##1lr##4##<t#<tJt#AA# 


STATION 

PHI 

K/( S*SQRT(PI A/Q) ) 

FORCE -METHOD rnn MFTunn 

********************************i,1 t 1cl,i<i e1t i t1c i'* i ,+**+^~~“*^"™ 

1 

0.000 

0.11946E+01 

0.10312E+01 

2 

11.250 

0.11 966E+01 

0. 11288E+01 

3 

22.500 

0.13205E+01 

0. 13348E+01 

4 

33.750 

0.14756E+01 

0. 15066E+01 

5 

45.000 

0.16298E+01 

0.16203E+01 

6 

56.250 

0.17546E+01 

0.17679E+01 

7 

67.500 

0. 18620E+01 

0. 18543E+01 

8 

78.750 

0.19784E+01 

0. 19386E+01 

9 

90.000 

0.21049E+01 

0.21615E+01 


ALL ELEMENTS SATISFY EQUILIBRIUM 


STOP (called by $MAIN ) 

CPr 45.308s, Wallclock: 116.253s, 9.7% of 4-CPU Machi 

HMM mem: 7733236, HWM stack: 310499, Stack overflows: 0 


ne 


Table 10: Output file outem28 for Example 6. 


m*********M*******************M* M **** M *** mwtMHm 

EMBEDDED CRACK-REMOTE TENSION A/C-0.2 ,A/T=0 8 
****************************************^; 4 ;^;^^^^ 


DESCRIPTION OF THE MODEL 


OUTPUT OPTION 
YOUNG S MODULUS 
POISSION S RATIO 

NUMBER OF NODES IN THE MODEL = 
NUMBER OF ELEMENTS IN THE MODEL- 

SHORT 

0.300000E+08 

0.300 

2464 

1856 

NODE 

NODAL COORDINATES 

X-COORD Y-COORD Z-COORD 

1 

5.00000 0.00000 

0.00000 
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2 

3 

4 


5.00266 

5.00245 

5.00187 


0.00000 

0.00510 

0.00930 


0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN- 0 


********************************************************************** 

********************************************************************** 


*********************************** 
equilibrium 
SUM OF THE X FORCE- 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


*********************************** 

CHECKS 

0.8550499E-09 
-0.4030978E-08 
-0. 1866076E-09 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 

X-COMPONENTS: FORCE- ° p AREA- 

v rnMPDNFNTS* FORCE* 0.62496ZE+OZ AKtH 

Z- COMPONENTS: FORCE- 0.000000E+00 ARE ^“ 


0.000000E+00 

0.624962E+02 

0.000000E+00 


nominal stresses 

NOMINAL STRESS IN THE X- DIRECTION - O.OOOOOOOE+OO 

NOMINAL STRESS IN THE Y- DIRECTION - 
NOMINAL STRESS IN THE Z- DIRECTION - 


0. 1000000E+01 
O.OOOOOOOE+OO 


************************************************************ 
STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************ 1 ********************************************* 


FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE-K 

K/(S*SQRT(PI*A/Q) ) 

1 

0.000 

0.9477525E+00 

0.5615390E+00 

2 

11.250 

0. 1026610E+01 

0.6082616E+00 

3 

22.500 

0. 1222969E+01 

0.7246035E+00 

4 

33.750 

0. 1440403E+01 

0.8534321E+00 

5 

45.000 

0. 1682060E+01 

0.9966126E+00 

6 

56.250 

0. 1921766E+01 

0.1138638E+01 

7 

67.500 

0.2169597E+01 

0. 1285476E+01 

8 

78.750 

0.2345304E+01 

0.1389582E+01 
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9 

90.000 

0.2409995E+01 

0. 1427911 E+01 


FROM THE CRACK OPENING DISPLACEMENT 

METHOD 

1 

0.000 

0. 76046 10E+00 

0.4505696E+00 

2 

11.250 

0.9681786E+00 

0.5736414E+00 

3 

22.500 

0.1228963E+01 

0.7281549E+00 

4 

33.750 

0.1465673E+01 

0.8684046E+00 

5 

45.000 

0. 1658270E+01 

0. 98251 74E+00 

6 

56.250 

0. 1928974E+01 

0. 1142908E+01 

7 

67.500 

0.2174445E+01 

0. 1288349E+01 

8 

78.750 

0.2351977E+01 

0. 1393536E+01 

9 

90.000 

0.2418940E+01 

0. 1433211E+01 


*************************ici'i'i' i ' i ' i ' i ' i ' i ' i ' ic1 ' it 


************************************************************ 


STATION PHI 

**************************** 


K/( S*SQRT(PI A/Q) ) 
FORCE-METHOD COD METHOD 

******************************** 


1 

0.000 

0.56154E+00 

0.45057E+00 

2 

11.250 

0.60826E+00 

0.57364E+00 

3 

22.500 

0.72460E+00 

0.72815E+00 

4 

33.750 

0.85343E+00 

0.86840E+00 

5 

45.000 

0. 99661 E+00 

0.98252E+00 

6 

56.250 

0. 11386E+01 

0. 11429E+01 

7 

67.500 

0. 12855E+01 

0. 12883E+01 

8 

78.750 

0. 13896E+01 

0. 13935E+01 

9 

90.000 

0.14279E+01 

0.14332E+01 


ALL ELEMENTS SATISFY EQUILIBRIUM 


STOP (called by SMAIN ) 

CP: 45.622s, Wallclock: 88.919s, 12 8% of 4-tpii Marhin« 

HWM mem: 7733236, HWM stack: 310499, SUck oJem oils: J 
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Table 11 : Output file occorl5 for Example 7 . 


^+++******************************************* 

*i2kL^***5I*i**i-**- 


DESCRIPTION of the model 


OUTPUT OPTION 

young s modulus 
KX S ntos°in the model 

Her OF ELEMENTS IN THE MODEL- 2260 


SHORT 

0.300000E+08 

0.300 

2863 


NODAL COORDINATES 


NODE 


X- COORD 


Y- COORD 


Z-COORD 


1 

2 

3 

4 


1.00000 

1.01320 

1.01220 

1.00930 


0.00000 

0.00000 

0.00520 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN- 0 

Iu8 “ot'tHe’A^E 0 - 0 381I9nE-09 

SUM OF THE Y FORCE- -0.6030213E 08 

SUM OF THE Z FORCE- -0.2428789E 08 

APPLIEOLOAD^d" THE SURFACE AREA COMPONENTS # 000000E+00 

X-COMPONENTS: FORCE- O.OOOOOOE+OO AKt^ 0<5 39980E+02 

"5SSS; force'- oioOOOOOEYOO AREA-_ O.OO0OOOE.OO_ 

nominal stress” in the x- OJ^CTION = o. iooooooi+oi 

«3£ S i N N THE Z- ffiS - 0 ; 0000000E.00___ 
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************************* 1 ,^^^ 


STATION PHI 

1 0.000 

2 11.250 

3 22.500 

4 33.750 

5 45.000 

6 56.250 

7 67.500 

8 78.750 

9 90.000 


FROM THE FORCE METHOD 

absolute-k 

0.2326628E+01 

0.2276628E+01 

0.2221275E+01 

0.221 8767E+01 

0.2265604E+01 

0.2362867E+01 

0.2530806E+01 

0.2691 461 E+01 

0.2359771E+01 


K/(S*SQRT(PI*a/qj ) 

0.2060499E+01 
0.2016218E+01 
0. 1967197E+01 
0. 1964975E+01 
0.2006455E+01 
0.2092592E+01 
0.2241322E+01 
0.2383601 E+01 
0.2089851 E+01 


1 

2 

3 

4 

5 

6 

7 

8 
9 


r*un iHt CRACK OPENING DISPLACEMENT METHOD 


0.000 
11.250 
22.500 

33.750 
45.000 
56.250 
67.500 

78.750 
90.000 


0.2387579E+01 
0.2243456E+01 
0.2194859E+01 
0.2191650E+01 
0.2237883E+01 
0.2333527E+01 
0.2482424E+01 
0.2787567E+01 
0.2379925E+01 


0.2114478E+01 
0. 1986840E+01 
0. 1943802E+01 
0.1 940961 E+01 
0. 1981904E+01 
0.2066608E+01 
0.2198474E+01 
0.2468713E+01 
0.2107699E+01 


NUMBER OF $IN6ULARITY*ELEMENTS*iN*THr*u******************** 
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******************************** *******************^^£*^£ 

K/( S*SQRT (PI A/Q) ) 

STATION PHI FORCE-METHOD COD METHOD 

A***************************************************^*^*^ 


1 

0.000 

0.20605E+01 

0.21145E+01 

2 

11.250 

0.20162E+01 

0. 19868E+01 

3 

22.500 

0. 19672E+01 

0. 19438E+01 

4 

33.750 

0.19650E+01 

0. 19410E+01 

5 

45.000 

0.20065E+01 

0. 19819E+01 

6 

56.250 

0.20926E+01 

0.20666E+01 

7 

67.500 

0.22413E+01 

0.21985E+01 

8 

78.750 

0.23836E+01 

0.24687E+01 

9 

90.000 

0.20899E+01 

0.21077E+01 


*********************************************************************£ 

LOADING NUMBER 2 

********************************************************************** 

********************************************************************** 
EQUILIBRIUM CHECKS 
SUM OF THE X FORCE- 0.7870824E-08 

SUM OF THE Y FORCE- -0. 1061700E-08 

SUM OF THE Z FORCE- -0.2249274E-06 


APPLIED LOAD AND THE SURFACE AREA COMPONENTS 


X- COMPONENTS: 

FORCE- 

O.OOOOOOE+OO 

AREA- 

O.OOOOOOE+OO 

Y- COMPONENTS: 

FORCE- 

-0.201877E-02 

AREA- 

0 . 539980E+02 

Z-COMPONENTS: 

FORCE- 

O.OOOOOOE+OO 

AREA- 

O.OOOOOOE+OO 


********************************************************************** 

NOMINAL STRESSES 
NOMINAL STRESS IN THE X- DIRECTION = O.OOOOOOOE+OO 

NOMINAL STRESS IN THE Y- DIRECTION - -0.3738595E-04 

NOMINAL STRESS IN THE Z- DIRECTION = O.OOOOOOOE+OO 


************************************************************ 

STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************************************************************ 


85 



FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE-K 

K/(S*SQRT(PI*A/Q) ) 

1 

0.000 

0. 1709044E+01 

0. 1513557E+01 

2 

11.250 

0. 1542926E+01 

0. 1366440E+01 

3 

22.500 

0. 1310887E+01 

0. 1160942E+01 

4 

33.750 

0.11 13279E+01 

0.9859379E+00 

5 

45.000 

0.9482749E+00 

0.8398074E+00 

6 

56.250 

0.8170198E+00 

0.7235657E+00 

7 

67.500 

0.7235576E+00 

0.6407942E+00 

8 

78.750 

0.6573300E+00 

0.5821419E+00 

9 

90.000 

0.5045022E+00 

0.4467952E+00 


FROM THE 

CRACK OPENING DISPLACEMENT METHOD 

1 

0.000 

0. 1814866E+01 

0. 1607274E+01 

2 

11.250 

0.1519294E+01 

0. 1345511E+01 

3 

22.500 

0. 1294912E+01 

0. 1146795E+01 

4 

33.750 

0. 1097889E+01 

0.9723082E+00 

5 

45.000 

0.9341272E+00 

0.8272779E+00 

6 

56.250 

0.8033078E+00 

0.7114222E+00 

7 

67.500 

0.7082324E+00 

0.6272219E+00 

8 

78.750 

0.6602082E+00 

0.5846909E+00 

9 

90.000 

0.5569405E+00 

0.4932353E+00 


************************************************************ 

NUMBER OF SINGULARITY ELEMENTS IN THE MODEL= 64 
************************************************************ 

************************************************************ 




K/( S*SQRT(PI A/Q) 

) 

STATION 

PHI 

FORCE -METHOD 

COD METHOD 

************************************************************ 

1 

0.000 

0.15136E+01 

0. 16073E+01 

2 

11.250 

0.13664E+01 

0. 13455E+01 
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3 

22.500 

0.11609E+01 

0. 11468E+01 

4 

33.750 

0 . 98594E+00 

0.97231E+00 

5 

45.000 

0. 83981 E+00 

0.82728E+00 

6 

56.250 

0.72357E+00 

0.71142E+00 

7 

67.500 

0.64079E+00 

0.62722E+00 

8 

78.750 

0.58214E+00 

0.58469E+00 

9 

90.000 

0.44680E+00 

0.49324E+00 


ALL ELEMENTS SATISFY EQUILIBRIUM 


Wallc^ock: 64.907s, 13.3% of 4-CPU Machine 

HUM mem: 7733236, HUM stack: 310499, Stack overflows: 0 


Table 12: Output file oscor!5 for Example 8. 


************************************************************ 
SURFACE CRACK AT A CIRCULAR HOLE A/C-l .A/T-0.5 , R/T=1.0 
************************************************************ 



DESCRIPTION OF THE MODEL 



OUTPUT OPTION 
YOUNG S MODULUS 
POISSION S RATIO 

NUMBER OF NODES IN THE MODEL * 
NUMBER OF ELEMENTS IN THE MODEL* 

SHORT 

0.300000E+08 

0.300 

2863 

2260 

NODE 

NODAL COORDINATES 

X- COORD Y-COORD Z- COORD 


1 

2 

3 

4 


1.00000 

1.01320 

1.01220 

1.00930 


0.00000 

0.00000 

0.00520 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


I ERR FROM SYMBN* 

********************************* 


************************************* 


LOADING NUMBER 


********************************* 


************************************* 
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equilibrium checks ********* 

SUM OF THE X FORCE* 0 .31 16301E-09 
SUM OF THE Y FORCE- -0.5856798E-08 
SUM OF THE Z FORCE- -0. 1453687E-09 


X- COMPONENTS 
Y-COMPONENTS 
Z-COMPONENTS 


i?c LIED L ?!°J ND the surface area’components 


FORCE- O.OOOOOOE+OO 
FORCE- 0 . 539980E+02 
FORCE- O.OOOOOOE+OO 


AREA- 

AREA- 

AREA- 


O.OOOOOOE+OO 

0.539980E+02 

O.OOOOOOE+OO 


**********************************^^^^^^^^ 

nominal stresses 

NOMINAL STRESS IN THE X- DIRECTION - o OOOonnnF^nn 

NOMINAL STRESS IN THE Y- DIRECTION - 0 ' 1 OOOOOOFtn? 

no«inal stress in the z- direction . O.OOOOOOOE+OO 


********************************^^^^^^^^^ 
STRESS INTENSITY FACTORS ARE AS FOLLOWS 
*******************************^ 


STATION 

PHI 

FROM THE FORCE METHOD 
ABSOLUTE-K 

K/(S*SQRT(PI*A/Q)) 

1 

0.000 

0. 1972444E+01 

0.1746828E+01 

2 

11.250 

0. 1981288E+01 

0.1 754661 E+01 

3 

22.500 

0.2008866E+01 

0. 1779084E+01 

4 

33.750 

0.2058871 E+01 

0. 1823369E+01 

5 

45.000 

0.2139874E+01 

0. 1895107E+01 

6 

56.250 

0.2259833E+01 

0.2001344E+01 

7 

67.500 

0 . 2442181 E+01 

0.2162834E+01 

8 

78.750 

0.2614491E+01 

0.2315435E+01 

9 

90.000 

0.2308667E+01 

0.2044592E+01 

1 

FROM THE CRACK OPENING DISPLACEMENT METHOD 
0.000 0 . 1951981 E+01 0.1728706E+01 

2 

11.250 

0. 1960570E+01 

0. 1736312E+01 

3 

22.500 

0.1987496E+01 

0.1 7601 58E+01 
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4 

33.750 

0.2035504E+01 

0.1802675E+01 

5 

45.000 

0 . 2 1 14568E+01 

0.1872695E+01 

6 

56.250 

0.2232217E+01 

0.1976887E+01 

7 

67.500 

0.2396266E+01 

0.2122171E+01 

8 

78.750 

0.2704875E+01 

0.2395481E+01 

9 

90.000 

0.2341106E+01 

0.2073321E+01 


I.i,************************* 1 '" 1 ' 1 ’—'' 

*5Sssl---_s! sasEsaas — sslsss? 


1 

2 

3 

4 

5 

6 

7 

8 
9 


0.000 

11.250 

22.500 

33.750 

45.000 

56.250 

67.500 

78.750 

90.000 


0.17468E+01 

0.17547E+01 

0.17791E+01 

0.18234E+01 

0.18951E+01 

0.20013E+01 

0.21628E+01 

0.23154E+01 

0.20446E+01 


0. 17287E+01 

0. 17363E+01 

0.17602E+01 

0.18027E+01 

0. 18727E+01 

0.19769E+01 

0.21222E+01 

0.23955E+01 

0.20733E+01 


ALL ELEMENTS SATISFY EQUILIBRIUM 


S° P « ( Q58 1 s ed Wall clock: 64.160s, 13.2% of 4-CPU Machine 

HWI menu ^7733236, HWM stack: 310499, Stack overflows: 0 
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Table 13: Output file osmcorl 5 for Example 9 ’ 


SURFACE CRACK AT A SEMI-CIRCULAR HOLE A/CM I/HT**.!** 


NODE 


description of the model 

OUTPUT OPTION y'imoT 

YOUNG S MODULUS , o°?Jnnnnr no 

POISSION S RATIO , 0.300000E+08 

NUMBER OF NODES IN THE MODEL I 2863 °’ 300 

NUMBER OF ELEMENTS IN THE HOTEL- 2260 

NODAL COORDINATES "" 

X-COORD Y-COORD 


Z- COORD 


1 

2 

3 

4 


1.00000 
1.01320 
1.01220 
1 . 00930 


0.00000 

0.00000 

0.00520 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


ICKK r KUrl SYMBN= 0 

’************‘****"* 7 ^“rN*^;R*****r**************‘*******‘***** 

rrrrrrrrr^rrrrrn 

SUM OF THE X FORCE* 0 1603649E 08 

fj* n°p thf 7 y ™ C t E - -O.ISI-OT 

?^.° F THE 1 F0RCEb -0.2318881E-09 


X- COMPONENTS 
Y-COMPONENTS 
Z- COMPONENTS 


APPUED LOAD AND THE SURFACE AREA COMPONENTS 
FORCE* 0. 000000 E+00 AREA* 

FORCE* 0.539980E+02 

FORCE- O.OOOOOOE+OO 


AREA- 

AREA- 


0.000000E+00 

0.539980E+02 

0.000000E+00 


***********************^^ a ^^ ## ^^ a ^^^ 

NOMINAL S T R F s 
NOMINAL STRESS IN THE X- DIRECTION * 

NOMINAL STRESS IN THE Y- D RECTI ON I 

NOMINAL STRESS IN THE Z- DIRECTION I 


********************* 
E S 

0.0000000E+00 

0.1000000E+01 

0.0000000E+00 
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** ********************************* 1 '+ 1 ' 1 ,* 1 ' 1 ' 1t **+*^^^^' k1 ' i ' i ' 


STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************************,Mr**************^****************** 


FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE-K 

K/(S*SQRT(PI*A/Q)) 

1 

0.000 

0.2132622E+01 

0. 1888684E+01 

2 

11.250 

0.2142476E+01 

0. 1897411E+01 

3 

22.500 

0.21731 18E+01 

0. 1924548E+01 

4 

33.750 

0.2228355E+01 

0. 1973466E+01 

5 

45.000 

0.231 70 19E+01 

0.2051989E+01 

6 

56.250 

0.2446855E+01 

0.2166974E+01 

7 

67.500 

0.2641224E+01 

0 . 23391 10E+01 

8 

78.750 

0.2817368E+01 

0.2495106E+01 

9 

90.000 

0.2474677E+01 

0.2191613E+01 


FROM THE 

CRACK OPENING DISPLACEMENT METHOD 

1 

0.000 

0.2110524E+01 

0. 1869114E+01 

2 

11.250 

0.2120098E+01 

0. 1877593E+01 

3 

22.500 

0.2150030E+01 

0. 1904100E+01 

4 

33.750 

0.2203130E+01 

0. 1951127E+01 

5 

45.000 

0.2289751E+01 

0.2027840E+01 

6 

56.250 

0.2417397E+01 

0.2140886E+01 

7 

67.500 

0.2592295E+01 

0.2295777E+01 

8 

78.750 

0.291 883 1 E+0 1 

0.2584963E+01 

9 

90.000 

0.2493143E+01 

0.2207967E+01 


a*********************************************^^^^^^^ 

NUMBER OF SINGULARITY ELEMENTS IN THE MODEL- 64 
***********************^***************** 11 * 11 **************** 
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************************************************************ 

K/( S*SQRT(PI A/Q) ) 

STATION PHI FORCE-METHOD COD METHOD 

************************************************************ 


1 

0.000 

0.18887E+01 

0. 18691E+01 

2 

11.250 

0.18974E+01 

0. 18776E+01 

3 

22.500 

0.19245E+01 

0 . 19041 E+01 

4 

33.750 

0.19735E+01 

0.19511E+01 

5 

45.000 

0. 20520 E+01 

0.20278E+01 

6 

56.250 

0.21670E+01 

0.21409E+01 

7 

67.500 

0.23391E+01 

0.22958E+01 

8 

78.750 

0.24951E+01 

0.25850E+01 

9 

90.000 

0.21916E+01 

0.22080E+01 


ALL ELEMENTS SATISFY EQUILIBRIUM 
STOP (called by SMAIN ) 

CP: 33.899s, Wallclock: 60.821s, 13.9% of 4-CPU Machine 

HUM mem: 7733236, HUM stack: 310499, Stack overflows: 0 


Table 14: Input file datl2d for Example 10. 


SURFACE CRACK -Prescibed displacements- A/C=1.0 A/T=0.2 

SHORT 

0. 30000 E+08 0.30000E+00 

2161 1664 

1 1.000000000 0.000000000 0.000000000 

2 1.013200000 0.000000000 0.000000000 

3 1.012200000 0.005200000 0.000000000 

4 1.009300000 0.009300000 0.000000000 

5 1.005200000 0.012200000 0.000000000 


2156 25.000000000 10.000000000 -5.000000000 

2157 25.000000000 15.000000000 -5.000000000 

2158 25.000000000 25.000000000 -5.000000000 

2159 25.000000000 45.000000000 -5.000000000 

2160 25.000000000 85.000000000 -5.000000000 

2161 25.000000000 125.000000000 -5.000000000 

1 210 1 2 211 210 1 3 212 1 

2 210 1 3 212 210 1 4 213 1 
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1658 2141 

1659 2142 

1660 2143 

1661 2144 

1662 2145 

1663 2146 

1664 2147 

1 

2 


2140 2154 

2141 2155 

2142 2156 

2143 2157 

2144 2158 

2145 2159 

2146 2160 


0 

0 


2155 

2156 

2157 

2158 

2159 

2160 
2161 

0 

0 


2085 2084 2098 2099 

2086 2085 2099 2100 

2087 2086 2100 2101 

2088 2087 2101 2102 

2089 2088 2102 2103 

2090 2089 2103 2104 

2091 2090 2104 2105 


0 

0 

0 

0 

0 

0 

0 


1877 

1 

0 

0 


1878 

1 

0 

0 


1879 

1 

0 

0 


1880 

1 

0 

0 


1881 

1 

0 

0 


2161 

0 

0 

1 


0 

1 

0 

0 

0 


REMOTE 

0 

0 

0 

0 

0 

0 

0. 

0000 

0.0000 0.0000 

203 

0 

1 

0 

O.OOOOOE+OO 

204 

0 

1 

0 

O.OOOOOE+OO 

205 

0 

1 

0 

O.OOOOOE+OO 

206 

0 

1 

0 

O.OOOOOE+OO 

207 

0 

1 

0 

O.OOOOOE+OO 

208 

0 

1 

0 

O.OOOOOE+OO 

209 

0 

1 

0 

O.OOOOOE+OO 

412 

0 

1 

0 

O.OOOOOE+OO 

413 

j i j 

0 

1 

0 

O.OOOOOE+OO 

414 

4 i r 

0 

1 

0 

O.OOOOOE+OO 

415 

0 

1 

0 

O.OOOOOE+OO 

416 

0 

1 

0 

O.OOOOOE+OO 

417 

0 

1 

0 

O.OOOOOE+OO 

418 

0 

1 

0 

O.OOOOOE+OO 

621 

0 

, 1 

0 

O.OOOOOE+OO 

622 

0 

1 

0 

O.OOOOOE+OO 

623 

0 

1 

0 

O.OOOOOE+OO 

624 

0 

1 

0 

O.OOOOOE+OO 

625 

0 

1 

0 

O.OOOOOE+OO 

626 

0 

1 

0 

O.OOOOOE+OO 

627 

0 

1 

0 

O.OOOOOE+OO 

830 

0 

1 

0 

O.OOOOOE+OO 

831 

AAA 

0 

1 

0 

O.OOOOOE+OO 

832 

0 

1 

0 

O.OOOOOE+OO 

833 

0 

1 

0 

O.OOOOOE+OO 

834 

AA p 

0 

1 

0 

O.OOOOOE+OO 

835 

0 

1 

0 

O.OOOOOE+OO 

836 

0 

1 

0 

O.OOOOOE+OO 

1039 

0 

1 

0 

O.OOOOOE+OO 

1040 

0 

1 

0 

O.OOOOOE+OO 

1041 

0 

1 

0 

O.OOOOOE+OO 

1042 

0 

1 

0 

O.OOOOOE+OO 


0. 10000 E- 05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 1 0000E - 05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 
0.10000E-05 
0. 10000E-05 
0. 10000E-05 
0. 10000E-05 


O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

0 . OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 
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1043 0 1 0 

1044 0 1 0 

1045 0 1 0 

1248 0 1 0 

1249 0 1 0 

1250 0 1 0 

1251 0 1 0 

1252 0 1 0 

1253 0 1 0 

1254 0 1 0 

1457 0 1 0 

1458 0 1 0 

1459 0 1 0 

1460 0 1 0 

1461 0 1 0 

1462 0 1 0 

1463 0 1 0 

1666 0 1 0 

1667 0 1 0 

1668 0 1 0 

1669 0 1 0 

1670 0 1 0 

1671 0 1 0 

1672 0 1 0 

1875 0 1 0 

1876 0 1 0 

1877 0 1 0 

1878 0 1 0 

1879 0 1 0 

1880 0 1 0 

1881 0 1 0 

1923 0 1 0 

1937 0 1 0 

1951 0 1 0 

1965 0 1 0 

1979 0 1 0 

1993 0 1 0 

2007 0 1 0 

2021 0 1 0 

2035 0 1 0 

2049 0 1 0 

2063 0 1 0 

2091 0 1 0 

2105 0 1 0 

2119 0 1 0 

2133 0 1 0 

2147 0 1 0 

2161 0 1 0 

0 0 0 0 

1 

1 2 

6 7 


0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

0.00000E+00 

O.OOOOOE+OO 

0.00000E+00 

0.00000E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

3 4 

8 9 


0. 10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

0. 10000E-05 

0. 10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0. 10000E-05 

0. 10000E-05 

0. 10000E-05 

0. 10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

0. 10000E-05 

0. 10000E-05 

0.10000E-05 

0.10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

0. 10000E-05 

0. 10000E-05 

0.10000E-05 

0. 10000E-05 

0.10000E-05 

O.OOOOOE+OO 

5 

10 


O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 


125.0000 25.0000 0.2000 5.0000 


1.0000 
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T^M e" 1 5I Output' f i 1 i 0«tdl2 for Example 10. 

******************************************** 

surfac e*c^ck*in*a*plate-displacehents^v : 

description of the model 


OUTPUT OPTION . S 0°300000E+08 

young s MODULUS _ 0.300 

NUHBER°OF S NODES IN THE MODEL * 2161 

HE OF ELEMENTS IN THE MOOEL- 1664 

nodal coordinates 

x- COORD Y -COORD 1 


1.00000 

1.01320 

1.01220 

1.00930 


0.00000 

0.00000 

0.00520 

0.00930 


I ERR FROM SYMBN= 


********************* 


****************** 
LOADING NUMBER 


**************** 


******************** 


******************** 


REACTION FORCES AT PRESCRIBED DISPLACEMENT NODES 


SUM OF THE X FORCE- 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


0.3833682 E- 04 
0.2999872E+02 
-0. 1479277E-09 


******************************^ 

SUM OF THE Y FORCE* -0-2««72E 09 

SUM OF THE Z FORCE- 0.S562037E 10 

rrr.^ 


Q-Z' 



STATION 

1 

2 

3 

4 

5 

6 

7 

8 
9 


FROM THE FORCE METHOD 


1 

2 

3 

4 

5 

6 

7 

8 
9 


PHI 

0.000 

11.250 

22.500 

33.750 

45.000 
56.250 
67.500 

78.750 

90.000 


ABSOLUTE-K 

K/(S*SQRT(PI*A/Q)) 

0.3113099E+00 

0.2757010E+00 

0.3026873E+00 

0. 268064 7E+00 

0.2916434E+00 

0.2582840E+00 

0.2850461 E+00 

0.2524414E+00 

0. 28101 11E+00 

0.2488679E+00 

0.2786616E+00 

0.2467871 E+00 

0.2773719E+00 

0.2456450E+00 

0.2767214E+00 

0.2450689E+00 

0.2765203E+00 

0.2448908E+00 


from the crack opening displacement method 


0.000 

11.250 

22.500 

33.750 

45.000 
56.250 
67.500 

78.750 

90.000 


0. 3202891 E+00 
0.2983115E+00 
0.2885413E+00 
0.2821246E+00 
0.2782588E+00 
0. 27601 32E+00 
0.2747916E+00 
0.2741661 E+00 
0.2739710E+00 


0.2836531 E+00 
0.2641894E+00 
0.2555368E+00 
0.2498540E+00 
0.2464304E+00 
0. 244441 7E+00 
0.2433598E+00 
0.2428059E+00 
0.2426330E+00 


STATION pur S*SQRT(PI A/Q) ) 

************.**.^^ ### ^^^^FORCE-METHOO^^^^cOD^METHOO 


1 

2 


0.000 

11.250 


0.27570E+00 

0.26806E+00 


0.28365E+00 

0.26419E+00 


96 



3 

22.500 

0.25828E+00 

0.25554E+00 

4 

33.750 

0.25244E+00 

0.24985E+00 

5 

45.000 

0.24887E+00 

0.24643E+00 

6 

56.250 

0.24679E+00 

0.24444E+00 

7 

67.500 

0. 24565E+00 

0.24336E+00 

8 

78.750 

0.24507E+00 

0. 24281 E+00 

9 

90.000 

0.24489E+00 

0.24263E+00 


ALL ELEMENTS SATISFY EQUILIBRIUM 
STOP (called by (MAIN ) 

CP: 28.635s, Wallclock: 55.219s, 13.0% of 4-CPU Machine 

HWM mem: 7733236, HWM stack: 310499, Stack overflows: 0 


Table 16: Output file outdxl2 for Example 11. 


♦♦★♦a*****************************************************^ 

SURFACE CRACK IN A PLATE-PRESCRIBED DISPLACEMENTS - U= -0.3E-7 on x*b 
A/C-1.0 A/T-0.2 / on x d 


DESCRIPTION OF THE MODEL 


OUTPUT OPTION 
YOUNG S MODULUS 
POISSION S RATIO 

NUMBER OF NODES IN THE MODEL = 
NUMBER OF ELEMENTS IN THE MODEL* 

SHORT 

0.300000E+08 

0.300 

2161 

1664 

NODE 

NODAL COORDINATES 

X-COORD Y-COORD Z-COORD 


1 

2 

3 

4 


1.00000 

1.01320 

1.01220 

1.00930 


0.00000 

0.00000 

0.00520 

0.00930 


0.00000 

0.00000 

0.00000 

0.00000 


HEIGHT OF THE MODEL = 125.00 
WIDTH OF THE MODEL = 25.00 
SURFACE LENGTH OF THE CRACK = 1.00 
DEPTH OF THE CRACK = 1.00 
THICKNESS OF THE PLATE * 5.00 
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$ JSfiS : i:S 

*KSi*L2LI5L*I5***55*lJ2t**********2;2!******************* 


MAXIMUM BANDWIDTH 

TOTAL CORE REQUIREMENT OF BIGK- 


1002 

4449930 


************************************************************ 
SUM OF THE FORCES BEFORE BOUNDARY CONDITIONS 

*E25**2S2I2S**2!!2IXI2!!****x****l**************************** 


SUM OF THE X- FORCES ARE 
SUM OF THE Y- FORCES ARE 
SUM OF THE Z- FORCES ARE 


0 . OOOOOOE+OO 
O.OOOOOOE+OO 
0. OOOOOOE+OO 


PROJECTED SURFACE AREAS IN 
SURFACE AREA X-COMPONENT 
SURFACE AREA Y- COMPONENT 
SURFACE AREA Z-COMPONENT 


EACH OF THE COORDINATE DIRECTIONS 
O.OOOOOOE+OO 
O.OOOOOOE+OO 
O.OOOOOOE+OO 


VOLUME OF THE SOLID MODELED = 0.156250E+05 


AT SUBPROGRAM SOLVE-B CPU TIME’ 
ACCUMULATED CPU- 0.453569E+01 


AT SUBPROGRAM SOLVE- 
ACCUMULATED CPU- 


E CPU TIME- 
0.276978E+02 


0.453569E+01 


0.231621E+02 


I ERR FROM SYMBN- 0 


******************************* 


*************************************** 


LOADING NUMBER 


********************************** 


************************************ 


REACTION FORCES AT PRESCRIBED DISPLACEMENT NODES 


SUM OF THE X FORCE- 
SUM OF THE Y FORCE- 
SUM OF THE Z FORCE- 


-0.2250000E+02 

-0.6472670E-08 

0.1818062E-10 


********************************************************************** 
EQUILIBRIUM checks 
SUM OF THE X FORCE- 

SUM OF THE Y FORCE- -0 - 2951272E- 10 

SUM OF THE Z FORCE- 


-0.8830405E-11 
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************************************************************ 

STRESS INTENSITY FACTORS ARE AS FOLLOWS 
************************************************************ 

FROM THE FORCE METHOD 


STATION 

PHI 

ABSOLUTE-K 

K/(S*SQRT(PI*A/Q) ) 

1 

0.000 

-0. 1037153E-04 

-0.9185194E-05 

2 

11.250 

-0.2597900E-05 

-0.2300741E-05 

3 

22.500 

0.9410255E-06 

0.8333872E-06 

4 

33.750 

0.5019746E-05 

0.4445567E-05 

5 

45.000 

0.9920509E-05 

0.8785760E-05 

6 

56.250 

0. 1478098E-04 

0. 1309027E-04 

7 

67.500 

0. 1890774E-04 

0. 1674499E-04 

8 

78.750 

0.2166702E-04 

0. 1918866E-04 

9 

90.000 

0.2263688E-04 

0.2004758E-04 


FROM THE 

CRACK OPENING DISPLACEMENT METHOD 

1 

0.000 

-0.9245535E-04 

-0.8187993E-04 

2 

11.250 

-0.7967842E-04 

-0.7056447E-04 

3 

22.500 

-0.6818843E-04 

-0.6038876E-04 

4 

33.750 

-0.4885792E-04 

-0.4326935E-04 

5 

45.000 

-0.2653945E-04 

-0.2350376E-04 

6 

56.250 

-0.4131448E-05 

-0.3658876E-05 

7 

67.500 

0. 1484371E-04 

0.1314583E-04 

8 

78.750 

0.2753143E-04 

0.2438228E-04 

9 

90.000 

0.3198674E-04 

0.2832797E-04 
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